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ABSTRACT 

Future space-based gravitational-wave detectors, such as LISA/SGO or a similar European 
mission (eLISA/NGO), will measure the masses and spins of massive black holes up to very 
high redshift, and in principle discriminate among different models for their evolution. Be- 
cause the masses and spins change as a result of both accretion from the interstellar medium 
and the black-hole mergers that are expected to naturally occur in the hierarchical forma- 
tion of galaxies, their evolution is inextricably entangled with that of their galactic hosts. 
On the one hand, the amount of gas present in galactic nuclei regulates the changes in the 
black-hole masses and spins through accretion, and affects the mutual orientation of the spins 
before mergers by exerting gravito-magnetic torques on them. On the other hand, massive 
black holes play a central role in galaxy formation because of the feedback exerted by AGN 
activity on the growth of structures. In this paper, we study the mass and spin evolution of 
massive black holes within a semianalytical galaxy-formation model that follows the evolu- 
tion of dark-matter halos along merger trees, as well as that of the baryonic components (hot 
gas, stellar and gaseous bulges, and stellar and gaseous galactic disks). This allows us to study 
the mass and spin evolution in a self-consistent way, by taking into account the effect of the 
gas present in galactic nuclei both during the accretion phases and during mergers. Also, we 
present predictions, as a function of redshift, for the fraction of gas-rich black-hole mergers 
- in which the spins prior to the merger are aligned due to the gravito-magnetic torques ex- 
erted by the circumbinary disk - as opposed to gas-poor mergers, in which the orientation of 
the spins before the merger is roughly isotropic. These predictions may be tested by LISA or 
similar spaced-based gravitational-wave detectors such as eLISA/NGO or SGO. 
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1 INTRODUCTION 

Massive black-hole (MBH) 
be the brightest sources of 
ture space-based detectors such as 
ferometer Space Ante n na (LISA) 



mergers are expected to 
gravitational waves for fu- 
the Laser Inter- 
jBender et al.1 1 19981; 



Danzmann Rudigerl 120031 ; IThe LISA International S cience Team! 
201 ll) or a similar m is sion le d by ESA (e LISA/ NGO, 
see lAmaro-Seoane et alj d20 1 2h : Ijennrich et al] d2012h ) or 
NASA (SGO. see IThe SGO Core Concept TearrJ feoilh ). These 
detectors are expected to be capable of observing tens or even 
hundreds of mer ger events during their lifetime, up to redshift 
z ~ 1 or larger dSesana. Volonteri. & HaardtH2007l : ISesana et al.1 
1201 ll ; lArun et alj|2009t) . and measure black-hole masses and 
spins with astonishing accur acy (~ 10~ 3 for the masse s 
and 10~ 2 for the spins, see iBerti. Buonanno. & Wilil J2005h : 
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lLang & Hughesl feOOd l200l |2008|) ; lLang. Hughes, & Cornishl 
d201D ). Also, they should be able to tell a binary of black holes 
with aligned spins from one with misaligned (and therefore 
precessing) spins by lookin g at the higher-order harmonics of 
the gravitational waveforms dLang & Hughesl l200d |2007|. 120081 : 
lLang. Hughes. & CornisrJl201lh . 

Massive black holes are also expected to play a crucial role 
in galaxy formation: in fact, in Active Galactic Nuclei (AGNs) 
accretion onto the central MBH is believed to power jets or disk 
winds capable of exerting a feedback on the growth of structures, 
by ejecting gas from the int erstellar medium (ISM) and from the 
intergalactic medium (IGM) dGranatq et alJl2004lLapi et alfeood 
ICroton et alj|200d iBower et alJl200dlHopkins et aljl2008F there- 
fore quenching star formation in low-redshift, massive galaxies. In- 
deed, this feedback mechanism is a central ingredient of our cur- 
rent understanding of galaxy formation, because it helps explain- 
ing why large dark-matter halos present low baryonic masses, re- 
sulting in a sharp cutoff at the high-mass end of the stellar mass 
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function that is not observ ed in the halo mass function dBell et all 
l2003al ; lBenson et al]|2003h . Also, it helps making sense of the so- 
called "anti-hierarchical" evolution (or "downsizing") of baryonic 
structures, i.e. the fact that the most massive galaxies are domi- 
nated by old stellar populations, while low-mass galaxies generally 
present young stellar populations and longer-lasti ng star-formation 
activity dCowie et alj|l99q ; iFontanot et alj |2009). thus suggesting 
that massive galaxies assemble at higher redshift than low-mass 
ones. Naively, this behavior may seem in contrast with the "bottom- 
up"formation of dark-matter halos, which assemble hierarchically 
by a series of mergers, but it can be reproduced, at least in its main 
features, by semianalytic al galaxy-formation models inc luding the 
effect of AGN feedback dScannapieco. Silk. & Bo uwens 20Q5[). 

This link between MBHs and the larger-scale galactic prop- 
erties works also in the opposite direction, because it is the 
amount of cold gas present in galactic nuclei that regulates 
accretion onto the MBHs, and therefore their mass and spin 
evolution. Also, the 100-pc scale circumbinary di sks that are 
thoug ht to form after gas-rich mergers of galaxies dMaver et al.l 
|2007|) exert torques on the spins of the MBHs, aligning them 
by the time the binary's separation has shrunk to ~ 0.1 
pc iBardeen & Pettersonl Il975l; Bogdanovic. Reynolds. & Milled 
120071 ; IPerego et alj 120091 ; IDotti et alj l2010al : iMaio et al.l I2OI2E 
and further spin alignment occurs from this separ ation to the 
merg e r as a result of Post-Newtonian resonances ( Sch nittmanl 
|2004 iKesden. Sperhake. & BertH l2010h . Not only do these ef- 
fects profoundly influence the spin and mass evolution - be- 
cause if the black-hole spins are alig ned before the merger, the 
spin of the final bl a ck hole is la rger dTichv & Marronettil 12001 
Lousto et alj l20l"ol; iRezzolla et 



2009: 



Barausse & Rczzoll 



Buonanno, Kidder, & Lehnei 



20081 ; iKesdenl l2008h and 
l2008l ; iReisswig et aiTl2009l ; 



its mass low er (Tichy & Marronetti 

lLoustoetal]|2010[ iKesdenl 120081) - but they also affect whether 
the MBH resulting from the merger remains in the galaxy or 
gets ejected. In fact, numerical-relativity simulations of black- 
hole binaries have shown that mergers can produce final black 
holes with large kick velocities relative to the center of mass 
of the initial binary configuration, due to anisotropic emission 
of gravitational waves. In particular, for equal-mass configura- 
tions with initial black-hole spins lying on the orbital plane 
of the binary, the kick velocity can be as large a s 2500-4000 
km/s (iCampanelli et al. 2007; Gonza lez et alj l2007h . and veloci- 
ties as high as 5 000 km/s may occur for configurations with off - 
equatorial spins dLousto &"z iochowerl l201 ll ; lLousto et al.l l2012h . 
Such velocities are larger than the typical galaxy escape veloci- 
ties, thus leading to the staggering conclusion that most galaxies 
may not host a MBH. This would be in stark contrast with observa- 
tions that most galaxies do host a MBH at low redshifts, and it may 
be bad news for hierarchical galaxy formation models, which as 
mentioned above heavily rely on the feedback from MBHs. How- 
ever, for binaries with aligned spins the kick velocity is consider- 
ably lower and typically not sufficient for ejecting the final black 
hole from the galaxy. Because, as mentioned above, circumbinary 
disks tend to align the black-hole spins prior to mergers, it is clear 
that tracking the evolution of gas in galaxies is crucial to correctly 
predict how many galaxies host a MBH, let alone the MBH mass 
and spin evolution. 

What makes galaxy formation a difficult problem is the huge 
range of scales that are involved, which go from the Gpc scale 
of the present cosmological horizon and the Mpc scale of typical 
z ~ galactic halos, through the 10-kpc scale of galactic disks 
and kpc scale of galactic spheroids, down to the 100-pc scale of 



circumbinary disks and pc scale of MBH accretion disks, and 
finally to the 10 -6 -10~ r pc scale at which MBH mergers take 
place. This huge dynamical range, together with the complex 
nature of the processes, often nonlinear and dissipative, that take 
place on small scales ("subgrid physics"), makes the problem very 
difficult to solve numerically in full generality. In fact, while the 
current paradigm of cosmological structure formation (the ACDM 
model) has enjoyed remarkable success in reproducing large-scale 
observati ons [such as the cosmic microwave background fluc- 
tuations jjarosik et alj|2011 ; Larson et alj 201 l b; the large scale 
cluster ing of galaxies dKomatsu et alj d2011 ); lEisenstein et alj 
(2005), and references therein); th e cosmic shear field measured 
through weak gravitational lensing dFu et al. (2008) and references 
therein); the small scale power spectrum of Lyman-alpha forest 
sources |jena et ai]|2005h ; the number density of galaxy clusters 
dHenrv et alj d2009f) and references therein)], a global understand- 
ing of galaxy formation can pres ently be attempted only by means 
of s e mi-analytica l mo de ls dK auffirmim. White, & Guiderdoni 
19931; IColeetalJ 1 19941 . |2000|; ISomerville & Primack! 1 19991; 
Somerville et alj 120081; ICroton et alj |2006|; IB ower et al. 
iBenson & Bowed |2010| ; iMonaco. Fontanot. & Ta ffoni 
or more ambitiously w ith h ydr odynamic al sim- 

1 Springel & Hemouistl 120031; ISpringell 2005 ; 

Di Matteo, Springel, & Hernquist 2005; Scannapieco et al. 2009; 
Dotti et alj|2007l; IDotti Colpi. & Haardj 120061; fpotti et alj Uoogl 
2010al ; IMaio et alj|2012l; lO'Shea et alj|2004l; iKravtsov & Gnedinl 
20051 iTassis et al. 



ulations 



20031; ITassis. Kravtsov. & Gnedird [2OO: 



Dubois & Teyssier 2008; Blecha et alj|201ll; iGuedes et aljhoill 



Siiacki, Springel, & HaehnehfeOl lLl2009l ; lTevssier et alj201lb" 

In this paper we use a semianalytical model for the forma- 
tion of galaxies in a ACDM universe to study the evolution of the 
spins and masses of MBHs in a self-consistent way, taking into ac- 
count both the feedback of the MBHs on the growth of structures 
and the influence of the galactic nuclear gas on the MBH accre- 
tion history and on the spin-alignment prior to mergers (which in 
turns affects the spin evolution, as well as the kick velocities of 
MBHs and thus their possible ejection from galaxies). For the pur- 
pose of this investigation, we adopt the widely accepted scenario 
in which the IGM collapses into disk structures, which give rise 
to bulges ("spheroids") when disrupted by major (i.e. comparable- 
mass) galactic mergers, or when they become self-gravitating and 
develop bar-instabilities. To describe the evolution of dark-matter 
ha los, we use a full extended Pres s-Schechter merger tree based 
on lParkinson. Cole. & Helfvl d2008l) . and then evolve the baryonic 
components along the merger-tree branches, employing analytical 
prescriptions at the "nodes" of the tree to mimic the effects of the 
mergers, and taking into account environmental effects (tidal strip- 
ping, tidal evapora tion and dynamical friction) using the results of 
lTaffonietalJd2003h . 

The MBH are evolved within this model starting from 
two possible scenarios for their seeds, namely a light-seed 
(M sec d ~ 150M©) scenario - in whic h the MBHs form as 
remnants of Pop III stars at z ~ 20 (Ma dau & Reesl l200lh 
- and a heavy-seed (M see d ~ 1O 5 M0) scenario - in which 
MBHs form from the c ollapse of massive protogalactic disks 
15 ( Koushiappas, Bullock, & Dekel 2004 



at 10 < 



iBegelman Volonteri. & Reesl 120061 ; lLodato & Nataraianl 12006 
The MBHs are then evolved along th e "branches" of th e 
merger tre e , usin g the "reservoir" model of iGranato et alj d2004l) ; 
lLapi et ail (2006) to describe the gas present in the nuclear re- 
gion, and taking into account the feedback they exert on the 
growth of structure by means of accretion-powered jets. As for the 
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accret ion mechanism, we follow iDotti et all J2010bh : iMaio et all 
d2012h and assume that the nuclear cold gas accretes coher- 
ently (i.e., with a fixed angular momentum direction) onto the 
MBH in a gas-rich enviro nment, thus res u lting in a stea dy in- 
crease of the MBH spin dBardeenl Il970l : iThornd 1 19741) . In a 
gas-poor environment, due to the absence of a rotationally sup- 
ported structure, we assume that the MBH accretes chaotically 
(i.e. in lumps of material with essentially random orientations 
of the orbital angular m omentum), which resu lts on average in 
a decrease of its spin fang & Pringld 1 2006b . Similarly, prior 
to a MBH merger, we assume that the MBH spins are aligned 
due to the gravito-magnetic torques exerted by the circumbi- 
nary disk if the nuclear environment is gas rich ("wet merger"), 
whereas we assume that they are randoml y oriented in a gas- 
poor nuclear environment ("dry merger") (Bardeen & Petterson 
19751: iBogdanovic. Reynolds, & Miller 20071; IPerego et al J 120091 
Dotti et al]2010ahlMaio et alj|2012h . 

We stress this idealized picture will be at least partially mod- 
ified by processes such as star formation and feedback, which can 
inject energy and angular momentum into the circumnuclear disk 
and potentially break t he coherence of t he accretion flow. However, 
recent simulations by Mai o et al. I J2012h show that these effects are 
probably not sufficient to create strongly turbulent environments 
and completely randomize the orbits of the accreting gas in the 
case of quasars fed by galaxy mergers, and that in these systems 
the accretion flow retains, at least partially, its coherence. The situ- 
ation could be different for quasars fed by disk instabilities at high 
redshifts, where cold flows a re important and accre t ion may take 
place b y massive clumps (e.g. lBournaud et alld201ll) ; lDubois et al.l 
( 1201 II) ). However, these clumps are expected to infall from large- 
scale disks, and might retain at least partially the disk's angular mo- 
mentum down to the MBH. In fact, observ ational hints for at least 
part ial alignment come from the re sults of lLagos et al.l d201lh (see 
also lLagos^ Padilla. & Cor"3 (2009)), who studied the relative ori- 
entation between AGNs and their host galaxies in the Sloan Digital 
Sky Survey, finding that random orientations are ruled out at high 
confidence. 

The effect of MBH mergers is accounted for by using analyti- 
cal formulas reproducing the results of numerical-relativity simu- 
lations. In fact, while the latter are the only way to study black-hole 
mergers rigorously, these simulations are very time-expensive and 
are not a viable option to cover the whole parameter space of black- 
hole binaries. Recently, the interface between numerical and analyt- 
ical relativity has produced a number of approaches which employ 
a combination of post-Newtonian theory, general-relativistic 
perturbation theory and/or fits to numerical data to reproduce 
different aspects of black-hole binaries, such as the gravitational 
wa veforms (see for inst a nce t h e effective-one-body approach, 
e.g.|Buonanno & DamouJ dl999h : |Pamour. JaranowskL & S chafer 



(2008); IBarausse & Buon annd d20ld); |Pan et al.l koill) 

[San 



e.g. i Santamarfa et alj d2010l) ). 



or 
the 



the hybrid waveforms 

kick velocity teaker et all 120071 120081; Ivan Meter et all l20ld : 
ICampanelli et al.l2007l;lLousto & Zlochoweill2009l ,1201 ll) . the final 
mass 1 Tichv & Marronettil2008l ; lLousto et alj20ld: 



Reisswig et 



al. 



20091; lKesdenll2008l) and the final sp i n dTichy & Marronettill2008k 
Lousto et al j I2OIO1; iR ezzolla et al. 2008; IBarausse & Rezzollal 
2009] ; iBuonanno. Kidder. & Lehnerll2008h. For th is investigation, 



we use the formula of 



Barausse & Rezzolla (2009) for the spin, the 



formulas of lTichv & Marronettil d2008h;lReisswig et alj d2009l) for 



Mete r et alj d2010h for the kick 



the mass, and the formula of 
velocity. 

Because we track the evolution of baryonic structures, 



and in particular of the gas present in galactic nuclei, along 
the dark-matter merger trees, we can attempt to discrimi- 
nate between chaotic or coherent accretion onto the MBHs, 
and between aligned-spin and precessing-spi n MBH mergers. 
This i mproves upon previo u s models, e.g. Berti & Vo lonteril 
(2008) ; IVolonteri etalj d2005|) . iFanidakis et all d2010ll201ll) "and 
lLagos. Padilla. & Coral d2009f) . which considered either chaotic 
or coherent accretion, and either aligned or misaligned merg- 
ersQ We stress that LISA/SGO or eLISA/NGO will be able to 
test our model not only by measuring the MBH masses and 
spins as a function of redshift dSesana. Volonteri. & Haardtll2007l; 
ISesana et "aill20T l|), but for each merger event it should also 
be able to determine whether the MBH binary has precessing 
or aligned spins, by looking at the higher order harmonics of 
the gravitational waveforms dLang & Hughes! 120061. 12007L 120081; 
lLang. Hughes. & Cornishll201lh . 

This paper is organized as follows. In Sec. [2] we present our 
semianalytical galaxy-formation model in detail, focusing on dark 
matter in Sec. 12. 11 on the baryonic components in Sec. 12. 21 and on 
merger and environmental effects in Sec. 12. 31 In Sec.[3]we calibrate 
the free parameters of our model to reproduce existing observations 
at z — and at z > 0. In Sec. [4] we present our predictions for 
the character of MBH mergers (i.e. whether they involve aligned 
or precessing spins), while in Sec. [5] we study the evolution of the 
spins of MBHs as a function of redshift. In Sec. [6] we draw our 
conclusions and present plans for future work. 

Throughout this paper, we assume a flat ACDM cosmology 
with n DM = 0.227, fl b = 0.0456, Hp = 70.4km/(sMpc) and 
cr 8 = 0.809 dJarosik et al]|201 lhlLarson et alj20"TH ; lKomatsu et alj 

Eoiih . 



2 PHYSICAL MODEL 

2.1 The dark-matter merger trees 

For the dark-matter merger and accretion e volution, w e adop t 
the extended Press-Schechter formalism of Cole et alj ^2000), 



as modified by IParkinsonT Cole, & Hellvl [2008). This algorithm 
reproduces the statistical properties of the dark-matter merger 
trees produced with cosmological nume rical N-body simula- 
tions dSpringel et alj|2005l : ICole et alj[ 2008). More specifically, we 
start our merger trees at an initial redshift z = 20 in the case 
of a light-seed M BH scenario (in whic h MBHs form as remnants 
of Pop III stars ( Ma dau & Reed 1200 lh ). while in the case of a 
heavy-seed scenario, in which MBH s form from the collapse of 
massive protogalactic disks dKoushiappas. Bullock. & Dekell2004l ; 
iBegelman. Volonteri. & Ree s 2006; L odato & Nataraiarl2006t) . we 
start our merger trees at z — 15. In the light-seed scenario, 
we assume that a halo of total mass A/ V ir forming at 15 < 
z < 20 contains a black-hole seed with mass M se ed = 150Mq 
if M v i r > 1.1 x 10 7 /i _1 Mq. This corresponds to populat- 
ing halos coll apsing from the larg e -cr peaks of the primordial 
density field dMadau & Reed 1200 ll ; IVolonteri. Haardt. & Madaul 



1 Indeed, iBerti & Volonteril J2008t) ; IVolonteri etaU (2005) did not de- 
scribe the baryonic components and therefore did not try to distinguish 
between gas-poor and gas- rich nuclear environments. On the contrary, 
IFanidakis et"al] d201(ll2oTTI) and lLagos. Padilla. & Cor3 420091) modelled 
the baryonic components in det ail, but did not try to i nfer w hether accretion 
is chaotic or coherent. Also, [Fanidakis et al. (2010, 2pTll) always consid- 
ered randomly oriented spins prior to mergers. 
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2003)0 In the heavy-seed scenario, instead, we place black-hole 
seeds with mass M scc d — 10 Mp in halos forming at 10 < 
z < 15 with virial mass ( Koushiap pas. Bullock. & Dekell |2004| ; 



ISesana. Volonteri. & Haardtll2007l) 



M vir > 10 7 



l + z 
18 



-3/2 



A 
O04 



-3/2 



Mp 



(1) 



where A is the halo spin parameter that we will introduce shortly. In 
both scenarios, because little is known about the spin of the seeds, 
we choose the spin parameter abh = cJbh/ (GM^) of each black- 
hole seed randomly from a uniform distribution — 1 < abh < 1. 
We stress, however, that the predictions of our model, and in par- 
ticular those regarding the spins of MBHs, are qualitatively inde- 
pendent of this assumption as long as ones looks at MBHs of mass 
M^h > 3M soc d . This is because a black hole loses memory of its 
initial spin when it accretes a mass comparable to its own (i.e., if 
a black hole of mass Mbh accretes coherently, which as we will 
show is the case at high re dshifts, it s spin becomes maximal after 
accreting a mass < 2Mbh dBardeerJll970T) ). 

Also, when it forms, a halo of total mass M v i r is assumed 
to contain unprocessed hot gas (see Sec. |2.2t with mass A/hot = 
/coii-Afvir, where the baryonic collapse fraction / co u is given by 



/coll (Mvir, Z) = 



fb 



(l + 0.26M / (z)/M vi r) 3 ' 



(2) 



with fb = Qb/^im ~ 0.16 (with f2 m = Qb + S2dm)- The fil- 
tering mass as a function of redshift, Mf(z), accounts for the 
effect of the ionizing extragalactic UV background produced by 
massive stars and quasars, which is ab le to partially reduce the 
baryonic content in low-mass systems dGnedinl fooO). In partic- 
ular, we calculate Mf(z) u s ing th e equations in Appendix B of 
iKravtsov! Gnedin, & Klvpin (2004), assuming ^overlap = 11 and 
Zrpinn = 10 dLarson et al J|20 1 ll) (z mn . T \ a p and 2 r eio» respectively 
correspond to the redshift at which multiple HII regions overlap, 
and to the redshift at which most of the medium is ionized). Af- 
ter the initial formation redshift, additional hot gas is brought in 
from the IGM by the dark matter that accretes onto the halo, and 
we therefore assume 



M inf = /coll M v ir 



(3) 



where the baryonic collapse fraction / coU is given again, as a func- 
tion of redshift, by Eq. {2). 

The resolution AA/ of the merger trees (i.e., the mass scale 
below which matter is assumed to accrete on an existing halo rather 
than give rise to a merger) is chosen to keep the computational time 
to acceptable levels while following the dark-matter halos (and the 
baryonic components within them) to very high redshifts. In par- 
ticular, we set AM = mm(10~ 3 Af , 10 10 M Q ) x (1 + z)~ 3 - 5 , 
Mo being the final mass of the halo at z — 0. The value of AM at 
z — is comparable or smaller than the resolutions typ ically used 
by semianalytical galaxy form ation models (see e.g. ICole et all 
(2000); So merville et alJd200 8)). and t he redshif t depen dence is in- 
troduced following Volonteri, Haardt, & Madau (2003) in order to 
track the merger tree to high redshifts. However, because our sim- 
ulations are computationally rather expensive (especially for large 
virial masses), in order to further cap the computational time while 
ensuring a range of masses that is sufficient to allow both minor 



and major mergers at all redshifts, at each redshift step we also stop 
following the branches that have mass smaller than S x M max (z), 
where S = 0.01 and M ma x(z) is the mass of the most massive halo 
at that redshift. 

The virial radius r v ir of a halo is related to its mass M v i r 
by the standard relation Mvii = 47rr 3 ir p C rit A c /3, where p CI n 
is the critical density, and where the density contrast at virial- 
ization, A c , is calculated following iBrvan & Normal] dl998l) . The 
halo density is assumed to be d escribed by the fitting function of 
iNavarro, Frenk. & White! d 19971) (NFW) 



PNFwM = Ps ( — 

. * S 



r B 



(4) 



truncated at the virial radius of the halo. The scale radius r s of the 
NFW profile is related to the virial radius by the so-called concen- 
tration parameter c(z) = r v ir/?V Imposing that the total halo mass 
equal M v i r , one immediately obtains that the scale density p s is 
given by p s = M vit (z)/\4-Kr^f(c)], with 



(z)/[4nr 3 s f(c)] 
/(c) = ln(l + c) 



1+c 



(5) 



The concentration parameter has been s tudied by several 
author s dBullock et al l 1200 ll ; IWechsler et all |2002| ; IZhao et all 
l2003allbl: uMaccio et aljl2007rh and found to present a large scat- 
ter for a fixed halo mass, but to s cale generally wi th the 
halo's main-progenito r history. Following Bullo ck et aljfcooll) and 
IWechsler et alj 12002), we adopt here a scaling c(z) oc 1/(1 + z) 
for the main-progenitor history after the halo formation: 



C MPH ( 2 ) = maX 



CO 

l + z 



Cf 



(6) 



where the concentration at formation is Cf = 4.1 dWechsler et al 
2002) and the concentration at z = is set following |Mjtctio_et_al 
J2007I) : 



log 10 c = 1.071 - 0.098 



lop 



Myir,0 

Mr, 



12 



(7) 



The concentration of halos not belonging to the main-progenitor 
history, however, is not expected to be given by Eq. (|6j. In- 
deed, smal l er halo s are expected to be more concentrated (cf. 
IZhao et all (2003a b), as well as the environmental effects de- 
scribed by Bull ock et al.l d200ll) ). To account for this effect, we 
adopt the high-concentration limit of the expressions of lZhao et al.l 
d2003al) . which give c oc M^T 1)/(3a) (with a = 0.48) at a fixed 
redshift. Combining this scaling with Eq. (|6j we obtain the expres- 
sion for the concentration of a halo of mass M v i r (z) at redshift z: 



c(z,M vil (z)) 



CQ / Myir(z) 



l + z \M MPH (2; 



(a-l)/(3e») 



(8) 

For simplicity, in this expression we assume that the main- 
progenitor history is given b y Mmph(z) = M v i r ,o exp(— Safz), 
with S = 2 and a f = cj /c dWechsler et all2002l) FI 

The angular momentum of each halo is determined by the 



halo's spin parameter, defined as A 



JvirS^M; 



5/2 



G~ 



where _E v ir and J v i r are the total energy and angular momentum 



2 In our cosmology, a mass of 1.1 X 10 7 h -1 Mq corresponds to the cos- 
mological Jeans mass collapsi ng at z = 20 from the 3.5 cr-peaks of the 
primordial density field iVolonteri. HaardC & Madau 20031). 



3 Applying Eq. (8) using the main-progenitor history extracted from the 
merger tree under consideration would be more difficult to implement in 
our code, and it is not clear that this procedure would be more accurate than 
the simple one that we use in this paper. 
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of the halo. We assign a spin parameter to a halo with no pro- 
genitors drawing from a log-normal distribution with median value 
A = 0.039 and standard deviation a =< y/(\n\ — In A) 2 >= 
0.53 dCole & Lacevlll996l : ICole et alJ fcoOO). We then assume that 
this spin parameter remains unchanged along the cosmic history 
of the halo, except if it experiences a merger with a second halo 
of comparable mass (i.e. if M v i r 2/M v i r i > 0.3), in which case 
we randomize the spin of the resulting halo by drawing from the 
same log-normal distribution. (We will discuss this in more detail 
in Sec. 12. 31). 



2.2 The baryonic matter 

2.2. 1 The hot gas phase 

We assume that the hot unprocessed gas phase is isothermal at the 
halo's virial temperature T v i r , and in hydrostatic equilibrium within 
the NFW profile, such that 



with 



Phot (r) = po exp 



27 



1 



ln(l + r/r s ) 



r/r s 



8-7T pm p Gp a r a 



(9) 



(10) 



27 k B T viT 

where m p is the proton mass, p is the mean molecular mass and po 
is calculated by normalizing to the total hot-gas mass at the redshift 
under consideration. The hot gas cools on a timescale t coo i given, 
at each point of its density distribution, by the standard expression 

3phot(r)fc s T vir 



TcoolM = 



2/xm p n e (r) 2 A(T vil .,Z)' 



(11) 



where n e (r) is the electron number density (which is proportional 
to Phot(V )), A(T, Z) is the cooling fu nction given by the tabulated 
results o f ISutherland & Do pita ( 1993), and where we assume that 
the hot gas is unprocessed and therefore has primordial metallicity 

z = io- 3 z Q . 

If the cooling time of the hot gas is shorter than its free-fall 
time, given approximately, at each radius, by 



3tt 



32Gp NFW (r) 



(12) 



(where p NFW (r) is the NFW average density within a radius r), 
then the hot-gas phase cools "fast" and undergoes gravitational col- 
lapse. If instead the cooling time is longer than the free-fall time, 
then the cooling is "slow" and the hot gas cools down through a 
sequence of quasi-hydrostatic equilibrium states. We therefore as- 
sume that the hot-gas phase is transferred into a cold-gas phase with 
rate 



M co ii(z) = 47r / — -f — ^-Ar , 

Jo t coll (r,z) 



(13) 



where £ co ii(r, z) = max(f coo i(r, z),t dyn (r, z)). 

However, on top of this "classical" picture employed al- 
read y by early semianalytical galaxy-form a tion mode l s (see 
e.g. iKauffmann. White. & Guiderdonil Jl993h ; ICole et al.l dl994l 
2000); ISomerville & Primackl 11 999)). we also include the effects 
of cold accretion flows foekel & Birnboiml2006l : lDekel et alj2009l ; 
ICattaneo et alj[200^) . which have been shown to be the predomi- 
nant mechanism leading to the formation of low-mass systems. In 
halos with mass lower than the critical mass 



M c = M s hock max[l, 10 



1.3(z 



where Af s h oc k = 2 x 1O 12 M0 and z c — 3.2, we assume 
that all the gas accreted from the IGM is not shock heated 
to the halo's virial temperatur e , but streams in on a dynamical 
time toekel & Birnboiml 120061 : IBekel et al.l 120091 : ICattaneo etail 
2006), thus enhancing star formation at high redshifts relative to 
the scenario where the accreting gas is shock heated. From the 
point of view of our model, this is equivalent to assuming directly 
tcoil = idyn in d!3t . for halos with M v j r < M c . Finally, in or - 
der to mimic the effe ct of ram pressure (Book & Benson 2010) 
and c lumpy accretion (Dekel & Birnbo im 2008; B irnboim & Dekell 
2011), which are expected to quench the cooling of the hot gas at 
low redshifts in large halos, we set t co \\ to the Hubble time in halos 
with M v i r (z) > 1O 13 M and z < 2. 



2.2.2 Density profile of the baryonic structures 

For the growth of the baryonic structures, we here adopt the widely 
accepted scenario in which the hot-gas collapse gives rise to a disk 
of cold gas with mass Ma >gas , where star formation may occur and 
result in the formation of a stellar disk with mass Md, stars- These 
disks may be disrupted by major galactic mergers and bar instabil- 
ities, thus forming a gaseous bulge with mass M(, lgas and a stel- 
lar bulge with mass M(, )S tars- Stellar formation takes place in the 
gaseous bulge as well, further contributing to stars- Also, as we 
will explain below, we allow part of the gaseous bulge to flow into a 
disk-like reservoir with mass M r0 s, which feeds the central MBHs 
during the accretion phase, and which forms the circumbinary disk 
expected to surround black-hole binaries after galactic mergers. 

Assuming a dissipationless collapse of the hot gas into the 
gaseous disk, the angular momentum conservation allows one to 
relate the halo's virial radius and spin parameter to the radius of the 
disk. In particular, adopting an exponential surface-density profile 
for the gaseous disk, 



r /r 



■(») 



E d (r, z) = T, (z)e 
the scale radius is given bv iMo. Mao. & Whitdfl998l) 

A jd r v 



with 



3 W 



f 

Jo 2 



y/2 m d ^% 



-f r (X,c,m d ,j d ) . 



l/(l + c) 2 -21n(l + c)/(l + c) 



[c/(l + c)-ln(l + c)] 2 



and 



f r (X,c,m d ,jd) 



U — rdu 

VclTvir) 



(15) 



(16) 



(17) 



(18) 



(14) 



where V c (r) is the velocity profile of the composite system (bulges, 
reservoir, disks, hot gas and dark matter) and where m d and j d 
are the ratios between the total mass and angular momentum of 
the disk and those of the halo. More specific ally, we take m d = 
(M^ rs + M| as )/A/vir and assume j d = m d IMo. Mao. & White! 
Il998h . 

In the calculation of the disk's scale radius, we also ac- 
count for the adiabatic halo response, which affects the velocity 
profile V c entering Eq . d 1 8b . using the standard prescription of 
Blumenthal et al. (1986). In particular, the angular momentum con- 
servation implies that the halo contraction following the baryonic 
collapse is described by 

M i {r i )r i = M } {r f )rf, (19) 

where r» and r/ are the initial and final radius of the shell under 
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consideration; Mi{r%) is the mass of the dark matter and hot gas 
contained in a radius n, calculated with the initial mass distribu- 
tion (which is assumed to be given by the NFW density profile for 
both the dark matter and the hot gas); and Af/(r/) is the mass of 
the composite system (dark matter, hot gas, disks, bulges and reser- 
voir) contained in the radius 77. Moreover, the mass conservation 
ensures 

M f (r f ) = 

M d {rf) + M b (r f ) + M res (r/) + M ho t{r f ) + M DM (r/) = 
M d (r f ) + M b (r f ) + M ICS {r f ) + (1 - / ga t)M j (r j ) , (20) 

where Mx (r) is the mass of component X within a radius r, and 
where / ga j = M ga ,i/M vir (with Af ga i = M d + M b + M ICS , M d 
and Mi being the total, i.e. gaseous and stellar, disk and bulge 
masses). Assuming spherical collapse without shell crossing, we 
can adopt t he ansatz rf = Frj, w here V = const is the contrac- 
tion factor (Blumentha l et al J 19861) . Eqs. dl9l > and j20\ can then be 
solved numerically for T, which then allows one to include the ef- 
fect of the adiabatic contraction following the baryon collapse in 
the calculation of the velocity profile V c of the composite system. 

For the stellar disk, we assume an exponential surface- 
den sity profile with scale radiu s r;f ara = r^ as /2 follow- 
ing ISomerville et alj d2008h a nd IZavala et alj J2003h (see also 
Dutto n & van den Boscbl j2009l) . who show that the gaseous disk 
is theoretically expected to have a larger scale radius than the stel- 
lar disk). As for the gaseous and stellar bulges, we assum e that they 
are described by the Hernquist profile lHerna uis3ll99d) 



Pb(r) 



Ml 



rb 



2tt r(r + r b ) a 



* = stars, gas , 



(21) 



where the scale radius r b (which we assume to be the same for 
the stellar and gaseous compone nts) is related to the half-light ra- 
dius Reg by r b = i? cff /1.8153 faemquisd[l99(ih . Using the fit of 
IShen et alj ( 120031) 1 



-5.54 + 0.56 log 



f °rl°gio(^) >10.3 



1 for R c 



log 10 (i? cff ) = < 

-1.21 + 0.14 log 10 (^) 
we can express the scale radius in terms of the total bulge mass 

Mb = M(, iSt ars + Mb, gas . 

For simplicitjQ we also assume that the reservoir can be de- 
scribed by an exponential surface-density profile with scale radius 
r res proportional to the influence radius of the central MBH, i.e. 

r rcs = aGM bh /VX, with a « 100. 

2.2.3 Star formation, supernova feedback and disk instabilities 

Assuming that star formation in galactic disks only happens in- 
side dense molecular clouds , which are well t raced by the HCN 
luminosity l lWu et alj 120051 ; iGao & Solomon] |2004|) , we follow 

4 The reservoir's geometry is needed for instance to calculate the velocity 
V c of the composite system, entering e.g. in Eq. 1181 . in the calculation of 
the adiabatic halo contraction factor T [Eqs. H9\ and )20H , and in that the 
total gravitational potential <f> appearing in Eqs. i27\ and )30t . However, the 
specific choice of the reservoir's geometry does not impact our final results 
significantly, because of its small size relative to the other components. 



iBlitz & Rosolowskvl d200fj|) ; button & van den Boscbl j2009h and 
take the disk star-formation rate (SFR) to be proportional to the 
molecular-cloud surface density as traced by HCN, £ mo i,HCN: 



E s f r = EsfSmol^HCN , 



(22) 



where e s f = 13Gyr _1 . From this expression one easily obtains the 
total SFR by integrating over the surface of the gaseous disk. 

More specifically, we write £ mo i,HCN as the product of the 
HCN fraction -Rhcn = S mo i : HCN/S mo i with the total molecular 
surface density of the disk, E mo i, which in turn we write as E mo i = 
/moiEd, gas > /moi being the molecul ar fraction of the disk ' s gas. For 
Ruck we use the fitting relation of Blitz & Roso lowskvl d2006h . 

i?HCN = 0.1 x (1 + E mol /(2OOAf pc- 2 )) - 4 , (23) 

and we relate f mol = flmoi/Cftmoi + 1) (with i? mo i = 
E mo i/E a to m ) to the mid-plane press ure P mp of the disk, follow- 
ing again lBlitz & Rosolowskvl feOOfjh : 



-Rmol — 



4.3 x 10 4 



(24) 



where the pressure and the Boltzmann constant are in CGS units. 
For the mid-plane pressure of the disk, we assume llElmegreenI 
1 19891; IBlitz & RosolowskyfeOOril ; button & van den Boschll2009l) 



-Pmp = gGEj, (E 9 + (cr fl /cr s )E s ) , 



(25) 



where we use o g /o s = 0.1 (button & van den~B osch 2009). 

We stress that in high-mass (and thus high-density) galaxies, 
where the molecular fraction / mo i ~ 1, this star-formation pre- 
scription reduce s to the standa rd Schmidt-Kennicutt star forma- 
tion power law |^nMcut3[l99i), E afr = e s f[E dig as/(M pc~ 2 )] n 
with n — 1.4 and e s f = 1.6 x 10~ 4 Mokpc~ 2 yr -1 , whereas 
in low-density systems the exponent n approaches 2.84. As a re- 
sult, the star formation law that we adopt suppresses star forma- 
tion in low-mass galaxies, in accordance with observation s (see 
button & van den Boscbl d2009l) : IBlitz & Rosolowskvl d2006h for a 
detailed discussion). 

The feedback from supernova events is expected to eject cold 
gas from the disk. To do so, the energy released by supernova ex- 
plosions in the disk must be sufficient to unbind the cold gas. We 
therefore compare, at each radius, the amount of energy released 
by these explosions with the binding energy, and write the total 
amount of cold gas ejected from the system as 



where 



Ml N (z) 



EsN(r, z) 



2tt 



rE S N(V, z)dr , 



£SN-EsNT?SNE 5 f r (r, z) 

(f>(r,z) 



(26) 



(27) 



Here, tj>(r, z) is the binding energy per unit mass of the composite 
system (bulges, disks, reservoir, hot gas and dark matter), tjsn is 
the number of Type II supernovae expected per solar mass of stars 
forme(Jf], Es~n = 10 44 J is the kinetic energy released per super- 
nova event, and esn is a parameter ranging from to 1 and regulat- 
ing the efficiency with which the supernova energy is transferred to 
the cold gas. We stress that the supernova feedback is most effective 



5 Following iRomano et a?] <2005l) , we adopt a IChabrieJ fe003l) initial 
mass function (IMF) between O.OOIMq and 1 Mq, and a power law 



cL/V/ dm* oc m 4 



therefore steeper than the standard Salpeter IMF, 



between 1 M and 100 Mq. This IMF gives rj S N = 5 X 10~ 3 M; 
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in low-mass systems, which present shallower potential wells from 
which the cold gas can easily escape due to supernova explosions. 

Also, disks are known to develop bar instabili- 
ties when they become self-gravitating, thus getting 
disrupted and transferring material to the bulge com- 
ponents dChristodoulou. Shlosman. & Tohlinel 1 19951 ; 

Efstat hTou. Lake. & Negropontdl 19821) . We assume that a stellar or 
gaseous disk is stable if, respectively, 



Vc{2.2r d ) 
{GM* d /r* d y/* 



> Qc 



* = stars , gas , 



(28) 



where o fffT 1.1 and a*Z 

0.9 dchristodoulou, Shlosman, & Tohlinel 1 19951 ; 

lEfstathiou. Lake. & Negroponte 19821) . If the disk becomes 
unstable, we assume that it gets disrupted in a dynamical time 
and transfers its material (either stars or gas) to the corresponding 
bulge component. 

As for the gaseous bulge, we assume that star formation takes 
place on a dynamical timescale. More specifically, we assume that 
the SFR per spherical shell in the gaseous bulge can be calculated 
as 



dtpb, ,s . 2 pb 
— — (r, t) — ATir 
dr 



(29) 



where i gas (r) = y / 3n/(32Gpb, sas ) is the local dynamical time for 
the bulge gas. Eq. l |29t can then be integrated over all radii to give 
the total SFR in the bulge. As in the disk case, we assume that the 
star formation ejects cold gas as a result of supernova explosions, 
at a rate 



M, 



b,gas 



(*) = 



£SN-EsNT?SNd^t(r, t)/di 
<t>(r,t) 



dr , 



(30) 



where the efficiency esN is assumed to be the same as for the disk. 
Again, this feedback mechanism is most effective in low-mass sys- 
tems. 

Finally we note that the SFR rates i22\ and J29t do not ac- 
count for the mass returned to the cold-gas phase by short-lived 
stars in the form of processed material. To include this effect, we 
use the instant-recycling approximation and assume that a fraction 
R — 0.29 of the mass is instantly returned into the cold-gas phase[f| 
In particular, this implies that the effective SFRs regulating the evo- 
lution of the disks and bulges are 



M b SFR (t) = (l- J R) / ^M)dr, 



(3D 



Mf R (t) 



[1-R) J 2± sb (r,t)nrdr . (32) 



2.2.4 The evolution of MBHs 



Star formation in the bulge i s believed to force, e .g. b 
radiation drag memural 1200 ll; iKawakatu & Umemural 
iKawakatu. Umemura. & Moril l2003h . part of the bulge's cold 
gas onto a low angular momentum circumnuclear reser- 
voir, which feeds the central MBH s during the accretion 
phases dHaiman. Ciotti. &~O striker 2004), and which may be iden- 
tified with the circumbinary disks expected to surround MBH bi- 
naries after a galactic merger. Because star formation in the bulge 



happens in violent bursts triggered by disk instabilities (see previ- 
ous section) or by galaxy mergers (as we will explain in the next 
section), we identify this accretion mode with the quasar mode of 
MBHs. 

In particular, we assume that the growth rate of the reservoir 
is pro portional to the bulge SFR dGranato et al.l 120041 ; lLapi et al.l 
2006) and is given by 



A reB ip b (t) . 



(33) 



where A ICS is a free parameter of our model. The cold gas in this 
reservoir then becomes available to accrete onto the central MBH 
at a rate 

M res 



Mqso 



(34) 



where the timescale t accr is a free parameter regulating the infall 
of the reservoir gas into the (pc-scale) MBH accretion disk. The 
accretion of this gas then changes the MBH mass according to 



M bh ,QSO = M QSO (l - v( a bh)) ■ 



(35) 



where the spin-dependent efficiency »7(abh) measures the energy 
emitted in electromag netic radiation by t h e accretion disk. More 
specifically, we follow iDotti etalj d2010bl) ; lMaio etal.1 d2012h and 
assume that the accretion onto the MBH takes place coherentl y 
(i.e. with a fix ed angular momentum direction, cf. lBardeenl dl970h : 
iThornd l fl974l) ) in a gas-rich environment, where gravito-magnetic 
torques rapidly align the the disk's angular momentum with the 
spin of the MBH. □ (As mentioned in the introduction, this ideal- 
ized coherent accretion flow may be at least partially randomized 
by star formation and feedback effects in the circumnuclear disk, or 
by the formation of clumps in high-redshift disk galaxies.) In a gas- 
poor environment, due to the absence of a rotationall y supported 
structu re, we assume that the MBH accretes chaotically dDotti et aU 
l2010bl) (i.e. in lumps of material with ess entially random orienta- 
tions of the orbital angular momentum, cf. lKing & Prin gle (2006)). 
Identifying a gas-rich environment with one where Af res > Mbh, 
we assume that the efficiency is 



r?(a bh ) = 1 - £ ISCO (a b h) 



(36) 



for M rcs > Mbh, with E ISCO the specific energy at the anti- 
clockwise innermost stable circular orbit (ISCO) around a Kerr 
black hole with spin parameter abh ranging between —1 (ex- 
tremal spin pointing do wnwards) and 1 (extremal spin pointing up- 
wards) dBardeenfl 970). Identifying instead a gas-poor environment 
with one where A-/ rcs < Mbh, we assume 



V( a hh) = 1 



EL 



,(Obh) + -Eisco( _a bh) 



(37) 



for M ros < Afbh, where we have assumed that accretion has equal 
probability of happening in the clockwise or anti-clockwise di- 
re ctions. (Thi s is a simplified version of the prescription derived 
bv lKing etai]d2005h .) 

Because of the energy and angular momentum carried by the 
gas accreting onto the MBH, in a gas-rich environment (A/ rcs > 
Mbh) the spin parameter abh increases steadily under coherent ac- 
cretion: 



«bh h QSO — [-^isco( a bh) — 2a b h-E ISCO ( a bh)] 



Mqso 
Mbh 



(38) 



6 To c alculate this retu rn rate, we follow again iRomano et alj 120051) and 
adopt a lChabrieil J2003l) IMF between O.OO1M and 1 M Q , and a power 
law diV/dm* oc m^ 27 between 1 M e and 100 M . 



7 Th is is the so-called Bardeen-Petterson effect dBardeen & Pettersonl 
1975), which also plays a fundamental role during gas-rich MBH mergers, 
as we will explain shortly. 
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where again E ISCQ and L lsco are the specific energy and angular 
momentum at the anti-clockwise ISCO around a Kerr black hole 
with spin parameter abh ranging between — 1 (extremal spin point- 
ing downwards) and 1 (extremal spin pointing upwards). In a gas- 
poor environment (Af rcs < Mbh), we again assume that accre- 
tion can happen clockwise or anticlockwise with equal probabili- 
ties, and the spin parameter decreases (on average) under chaotic 
accretion, following 



the radio-mode accretion and feedback are continuous and remains 
efficient up to z — 0. 

We therefore assume that when t coo \ > tg (i.e., when the 
hot gas undergoes quasi-hydrostatic cooling, cf. Sec. 12.2.1b the 
MBHs a ccrete directl y from the hot-gas phase at the Bondi accre- 
tion rate dBondill 19521) 



M bh , radio = 47rA s p hot (GM b h) 2 /w 



3 

s ) 



(44) 



• chaotic 
a bh,QSO 



,(a.bh) + L lsco (-a bh ) 



- a bh [£; isco (abh) + -E ISCO (-a bh )] 



1 Mqsq 

J M bh 



(39) 



We stress that we do not restrict the MBH accretion rate to val- 
ues lower than the Eddington rate, i.e. we allow super-Eddington 
mass accretion. However, follow ing the t heory of thin accretio n 
disks ( Shak ura & Sunvaevl dl973h : see also Poutanen et al. ((2003)), 
we assume the MBH's bolometric luminosity to be 



^bh.QSO 



M, 



QSOC , 



1+ln 



77(a b h)M Q soc 2 



(40) 



Also, thin-disk accretion is be lieved to produce jet outflo ws due to 
the Blandford-Znajeck effect (Blandford & Znaieldll977l) . and the 
jet power is parameterized by ( Meier 200 ll) 



Jjet,QSO 



r , 1n i2.7 -1 ( Ct N- - 1 0.9 / m V ' 

fiet x 10 ere s m Q - 

Jl 6 V0.01/ 9 V 0.1 , 



x (l + l.la bh + 0.29a bh ) 



(41) 



where a is the disk's viscosity parameter (for which we 
assume a = 0.1), m 9 = M b h/(1O 9 M ), rh = 
MQSo/(22m 9 M0yr _1 ), and where /j ct is a "fudge" fac- 
tor parameterizing the uncertainties affecting Eq. J41t (e.g., the 
Blandford-Znajeck jet-outflow power scales quadratically with the 
magnetic field, which is poorly known). These jets are expected to 
exert a feedback on the hot-gas phase and on the bulge gaseous 
component. More specifically, we assume th at they eject hot gas 
and bulge cold gas from the system with rates dGranato et al.l2004h 

2 £jct,QSO Mb 

3 



fcQBO 

b.gas 



(42) 



(43) 



cr 2 Af hot + M b , gas 

2 £jct,QSO Mhot 

3 a 2 M hot + M b , gas ' 

with a = 0.65Vvi r . 

In addition to the thin-disk quasar accretion mode, MBHs 
are expected to quiescently accrete matter directly from the hot- 
gas phase, when that is in quasi-hydrostatic equilibrium, through 
a thick advection-dominated accretion flow (ADAF). This is usu- 
ally dubbed "radio accretion mode" and does not contribute sig- 
nificantly to the mass evolution of MBHs, because the accretion 
rate is much smaller than for the quasar mode. However, the ra- 
dio mode is believed to play an important role in galaxy formation 
because it produces jet outflows that are much more powerful than 
thos e that would be produced by a thin disk with the sam e accretion 
rate dCroton et alj2006l : lBower et al.l200dlMeiej200lh . This com- 
pensates the smaller mass accretion rate of the radio mode relative 
to the quasar mode, and enhances its feedback. Moreover, while the 
quasar feedback is triggered by starbursts in the bulge, and there- 
fore intermittent and important mostly at relatively high redshifts, 



where phot is the density of the hot gas in the center of the galaxy; 
v s is the sound velocity in the hot gas, which is of the order of 
the virial velocity V V a', and As is a constant that depends on the 
adiabatic index of the gas, with As « 1. 12 for an isotherm al gas. 
The bolometric luminosity is then given ( Maha devanll 1 9971) by the 
ADAF luminosity 



dio = 1.3 x 10 J 



where m = Mbh, radio/ (22 trig Mq yr _1 ) and /3 i s related to vis- 
cosity parameter a = 0.1 by a w 0.55(1 — /3) (Fanida kis et al.l 
l201ll) . Because the radio-mode accretion happens through an 
ADAF and not through a thin disk, the rate of change of the spin pa- 
rameter is different from the quasar-mode case. More specifically, if 
we assume that the accretion is geometrically spherical, no angular 
momentum is transferred to the black hole and the spin parameter 
decreases due to the mass increase: 



^bh, radio 



-2dbh 



M b 



M b 



(46) 



Finally, ADAF accretion is expected to produce much more pow- 
erf ul jets than th e quasar mode, and the jet power is parameterized 
by dMeieill200ll) 



r 1 



, , ln 4B.i -1 / a \ - 1 / m \ 
/jetX 10 ergs [—) m 9 ^—jg 

(0.55/ 2 + 1.5/a bh + abh)) 



(47) 



where again mg = M bh /(1O 9 M ) and m = 
Mbh,radio/(22 mg M© yr _1 ), and wh ere / a nd q are dimen- 
sionless quantities, defined precisely in iMeiej d200ll) . regulating 
the actual angula r veloc i ty and azimuthal magnetic field of the 
system. Following |Meiei| [2001) we set / = 1 and g = 2.3, but we 
include a "fudge" factor /j e t [assumed to be the same as in Eq. J4 1 H 
to account for the uncertainties in Eq. ( |47| > (e.g. the uncertainties in 
the magnetic field, on which the jet power depends quadratically, 
and the higher-order terms in the black-hole spin, which are 
neglected in the standa r d Blandford-Znajeck calculation, cf. 
Tanab e & Nagatakil d2008l) ; iTchekhovskov. Naravan. & McKinnevl 
d2010h ). Like in the case of the quasar mode, we assume that 
the je ts remove hot gas a nd bulge cold gas from the system with 
rates dGranato et aT]|2004h 



Ml 



M b 



M, 



hot — 



2 -£*jet , radio 

3 a 2 Mhot + Mf,, g a S ' 

^ "jetj radio 

3 c 



M hot + M,, gas 



(48) 
(49) 



Finally, we stress that while the supernova feedback is most 
effective in low-mass systems, which presents shallow potential 
wells, the quasar and radio-mode feedbacks are most important in 
massive galaxies, where the quasar and radio activity is most pro- 
nounced. 



© 0000 RAS, MNRAS 000, 000-000 



The evolution of massive black holes and their spins in their galactic hosts 9 



2.3 Mergers and environmental effects 

The prescriptions outlined in the previous sections allow us to de- 
scribe the evolution of the baryonic structures along the branches of 
the dark-matter merger trees. Therefore, in order to produce "bary- 
onic merger trees" describing the complete evolution of the baryons 
along the cosmic history, we only need a prescription to describe 
what happens at the "nodes" of the dark-matter merger trees, i.e. 
when halos merge. 

First of all, we should note that the "nodes" of the dark-matter 
merger trees correspond to the instant at which the smaller halo 
(the "satellite") enters the bigger one (the "host"). After entering 
the host, the satellite halo survives as a bound substructure (a "sub- 
halo") within the host. However, because of the gravitational in- 
teraction with the particles of the host halo, known as dynamical 
friction, the satellite gradually loses energy and angular momen- 
tum and slowly sinks to the center of the host. It is only when the 
satellite halo reaches the center of the host that the subhalo finally 
dissolves and the baryonic structures of the two halos merge. 

Because our dark-matter merger trees are not necessarily bi- 
nary (i.e., their redshift step is such that a "node" may correspond 
to the merger of more than two halos), at each node we consider 
the biggest halo as the host, and calculate the dynamical-friction 
timescales for the remain ing halos (the satellites) using the fit - 
ting formula proposed bv lBovlan-Kolchin. Ma. & Ouataerj J2008t) . 
which accurately reproduces the merger timescales of extended ha- 
los predicted by N-body simulations: 



lite's center, at which the mean density of the satellite is compara- 
ble to the host's mean density at the satellite's position r: 



(M hoBt /M B 

at / 

Ki7 ln(l + Af host /M sat ) 



exp [c e] 



Rvii 



with A = 0.216, b = 1.3, c = 1.9 and d = 1. Here, e = L/L C (E) 
and r c (E) are respectively the "orbital circularity" (the ratio be- 
tween the satellite's angular momentum L and the angular mo- 
mentum L c of a circular orbit with the same energy E as the 
satellite) and the "circular radius" (the radius of a circular orbit 
with the same energy E as the satellite). These quantities describe 
the initial conditions of the infall (i.e. on what orbit the satel- 
lite enters the host). In particular, e may vary from to 1 and 
is drawn from a normal distr ibution centered on e ~ 0.5 and 
stand ard deviation a ~ 0.23 (Tormer] 1 19971 ; IKhochfar & Burkerj 
2006), while r c is derived from the periastron radius r p , which N- 
body simulations suggest should be correlated wit h the circularity e 
(more specifically, we assu me r v = R vir e 2 17 : cf. lTormenl ( fl99% ; 
IKhochfar & Burked I2006T) ). 

We notice that Eq. d50t already includes environmental effects 
such as the tidal stripping and tidal heating, which cause the mass 
of the satellite to decrease while sinking to the host's center. This 
is because Eq. ( 150b fits the results of N-body simulations, where 
these effects are automatically included. However, Eq. ( 150b does 
not include the effect of the continuous accretion of dark matter 
and hot gas onto the host halo along cosmic time, and we therefore 
correct for this effect at each redshift step of our merger tree by 
rescaling the remaining dynamical friction time to account for the 
host having grown in the meantime. 

Even though environmental effects are included in the dynam- 
ical friction timescale 1501) . as we have mentioned they also have the 
effect of reducing the satel lite's mass wh ile it sinks in the host. In 
particular, following iTaffoni et alj ( 120031) , we account for the tidal 
stripping (i.e. the tidal truncation of the satellite's density profile 
due to the average tidal force exerted by the host halo) by cutting 
the dark-matter and hot-gas density profiles of the satellite at the 
tidal radius -Rtidai, corresponding to the distance, from the satel- 



-NPW / D \ 
Psat (-KtidalJ 



-NPW/ n 

Phost (r) . 



(51) 



More specifically, because the tidal stripping is most effective when 
the satellite is at the periastron of its orbit, we assume r = r v in 
this equation and perform the cut when the satellite first reaches the 
periastron. Besides the tidal stripping, the satellite also loses mass 
due to the tidal heating, i.e. the evaporation induced by the rapidly 
varying tidal forces near the periastron. We assume that this effect 
causes both the dark matter as well as all the baryon components 
to lose mas s with characteristic timescale calculated as in the Ap- 
pendix B of Taffoni et al. ( 2003), and we assume that this mass loss 
starts when the satellites first reaches the periastron. 

When the satellite has sunk to the center of the host, a dy- 
namical friction time tdf after it first entered the host, the satel- 
lite halo finally loses its identity, and the baryonic structures of 
the host and satellite merge as well. When such a merger hap- 
pens, if the ratio MJgf /M^° st between the host and satellite halo 
masses is sufficiently large, we assume that the merger between 
the dark-matter structures perturbs the spin parameter of the re- 
sulting composite halo. More specifically, as already mentioned 



in Sec. [XT] if AC /Mv™ > 0.3 we assign the final compos- 
ite halo a new spin parameter A drawn from the same distribu- 
tion used for newly formed halos, i.e. a log-normal distribution 
with median value A = 0.039 and standard deviation a =< 
(50) ^/(lnA-ln A) 2 >= 0.53 dCole & Lacevlll996l : ICole etaihooch . 

while if M^/Mjg 8 * < 0.3 we leave the host's spin parameter 
unchanged. 

Also, further complications to this picture arise at the nodes 
where two or more halos, some or all of which containing their 
own subhalos, meet. In this case, if /M^° st > 0.3 for any 
one of the satellite halos, we recalculate all the dynamical fric- 
tion timescales using Eq. l |50| l, where we assume a new value for 
the circularity e, drawn from the same distribution used when the 
satellites first entered the host, i.e. a normal distri bution centered 
on e « 0.5 and standard deviation a w 0.23 jTormen|[l997l: 
IKhochfar & BurkertH2006l) . If instead M^'/M^' < 0.3 for all 
the satellite halos, we calculate the dynamical friction times of the 
satellites and their subhalos in the host using the procedure that we 
just described, but leave the dynamical friction times of the subha- 
los of the host unchanged. This scenario corresponds to the intuitive 
picture in which the incoming satellites manage to perturb and ran- 
domize the orbits of the host's subhalos only if they are sufficiently 
massive compared to the host. 

We also recall that numerical s imula- 
tions dWalker. Minos. & Hernouistl 1 19961 ; iNaab & Burkertl 120031: 
iBournaud. Jog. & CombesT 20051) suggest that mergers where 
the mass ratio between the total baryonic masses is larger than 
~ 0.25 — 0.3 ("major mergers") disrupt the galactic disks and 
give rise to a spheroidal component, while "minor mergers" (i.e., 
with mass ratio between the total baryonic masses smaller than 
~ 0.25 — 0.3) do not destroy the galactic disks, although they 
may perhaps drive the growth of the bulge by disk instabilities 
(see Sec. |2.2.3t . We implement this scenario in our model by 
defining a merger as "major" if the ratio of the baryonic masses 
is larger than 0.25 (Mb ary on, sat/M b aryon, host > 0.25 with 

-^baryon = Md, stars + Afd,gas + Afb, stars + -Mb, gas + Af r cs), 

otherwise we define the merger as minor. We then assume that in 
major mergers the disks are destroyed and their masses are added 
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to the corresponding bulge components 

M d , gas = 

M d , stars = (52) 

M(,, gas = Mj^gas + M^g as 
A^b.stara = M b< atara + M djStars 

M rcs = Af r h c ° s st + At£f , 

while we assume that minor mergers do not affect the galactic mor- 
phology, and therefore simply add the satellite's disk and bulge ma- 
terial to the disk and bulge component of the host galaxy: 

M d , gas = Aftg. + Mj^ 

(53) 

M b 

, stars — .stars .stars 

Mb, sas = Mb°g* B + M^ gas 

M ros = M r h c ° s st + M% . 

Another effect of mergers is to cause significant starburst events 
in the merging galaxies. This is automatically accounted for in 
our model, because the disruption of the galactic disks following 
a major merger channels gas into the spheroidal component, where 
star formation is very efficient, because it happens on a dynamical 
timescale [cf. Eq. J29H . 



2.3.1 Black-hole mergers 

When the baryonic structures of the host and satellite merge, a dy- 
namical friction time td{ after the satellite first entered the host, 
the central MBHs do not coalesce right away, but form a binary 
system. The binary then continues to harden t hrough "slingshot" 
interactions ( Sas law. Valtonen. & Aarsettill 19741) . in which stars in- 
tersecting the binary are ejected at velocities comparable to the bi- 
nary's orbital velocity, thus increasing the binding energy of the 
binary system. However, the binary will soon eject all the in- 
tersecting stars, thus making the slingshot hardening inefficient. 
This will cause the binary to stall at a separation of ~ 1 pc, un- 
less other mechanisms intervene to make it decay to a separa- 
tion of ~ 0.01 pc, where gravitational-wave emission becomes 
strong enough to drive the binary's evolution to the merger in 
a timescale shorter than the Hubble time. Since there is not, at 
present, a generally accepted scenario to overcome the stalling of 
the binary's evolution at ~ 1 pc, this bo ttleneck is generally re- 
ferred to as "the final-parsec problem" (Merritt & Milosavljevic 
l2005t) . It is generally thought, however, that the final-parsec prob- 
lem is somehow solved in nature, because uncoalesced binaries 
would result in slingshot ejection of MBHs when additional MBHs 
are brought in by successive mergers, thus resulting in off-center 
MBHs that seem rare or non-existent and in too much scatter in the 
M — a relation ( Haehnelt & Kau ffmann 2002). Also, some possi- 
ble mechanisms that would harden the binary until gravitational- 
wave emission becomes important, possibly solving the final- 
parsec problem, have been proposed. For instance, the presence of 
a gaseous acc retion disk would harden th e binary on the viscos- 
ity timescale iArmitage & Natara ian 2002); or the supply of stars 
available fo r slingshot interaction may be r eplenishe d by star-star 
encounters dMilosavlievic & Merrittl [20031 : mil [2002), or as a re- 
sult of the triaxial gravitational potential th at one naturally expects 
in galaxies forming from major mergers jlVIerritt & Poonl [20041 ; 
iBerczik et a0200d : lKhan. Just. & Merrittll201 ltlPreto et aU201 ll) . 



Because of all these uncertainties, and because the MBH coa- 
lescence timescale is in any case expected to be small compared to 
the dynamical friction time t&i and its uncertaintie o we make the 
simplifying assumption that the MBHs merge at the same time as 
the baryonic structures, i.e. a dynamical friction time tdi after the 
satellite first entered the host. This approximation is adopted, to our 
knowledge, by virtually all semianalytical galaxy-formation mod- 
els proposed so far, and has the advantage of greatly simplifying 
the implementation of our model0 

When the two MBHs merge, we can determine the mass 
and spin of the resulting MBH remnant if we make some as- 
sumptions on the relative orientation of the spins at large sep- 
arations. As mentioned in the introduction, in the last few 
years numerical relativity has reached a level of maturity suffi- 
cient for simulating black-hole binaries with non-aligned spins 
in a vast region of the space of parameters. Phenomeno- 
logical formulas iTichv & Marronettil 120081; iLousto et alj |2010|: 
Rezzolla etail [2008; Barausse & Rezzolla 2009; Reisswi getal] 



20091 : Buonanno . Kidder. & Lehnerl [2008; Kesden 2008) based 



on Post-Newtonian theory, general-relativistic perturbation theory, 
symmetry arguments, as well as fits to the numerical-relativity re- 
sults, are not only capable of reproducing to high accuracy the sim- 
ulation results for observables like the final mass, the final-spin 
magnitude and orientation, and the recoil velocity, but also allow 
one to make sensible predictions for these quantities in regions of 
the parameter space where simul ations are still to o time-expensive 
to ensure complete coverage (see lRezzollal j2009h for a rev i ew). 

Here, we use the formula of iBarausse & Rezzollal (2009), 
which predicts the MBH remnant' s spin magnitude an d orien - 
tation to very high accur a cy (se e IBarausse & Rezzollal d2009h : 
iKesden. Sperhake. & Bertil d2010h ). In particular, the final-spin 
magnitude is 



|ifin| = (ij^p [l a i| 2 + \o-2\ 2 q 4 + 2|a 2 ||ai|g 2 cosa + 



2 (|ai| cos/3 + |a 2 |g 2 cos7) \l\q + \l\ 2 q' 



1/2 



(54) 



with 



S4 



(1 + g 2 ) 2 
s 5 f + t + 2 
TTq 2 



i| 2 + a,2| 2 g 4 + 2|ai||a2|<j 2 cosa) + 



|ai| cos/3 + \a,2\q COS7) 



(55) 



Here, q — Mbh,2/Afbh,i is the mass ratio; v = q/(l + q) 2 is the 
symmetric mass ratio; |ai| and ei2 1 are the initial spin magnitudes; 
a, P (7) are the angles (at large separation) respectively between 
the two spins and between spin 1 (spin 2) and the direction of or- 
bital angular momentum, L; and S4 = —0.1229 ± 0.0075, S5 = 
0.4537± 0.1463, to = -2.8904 ±0.0359, t 3 = 2.5763±0.4833 
and t 2 = -3.5171 ± 0.12080 



8 For instance, for a satellite in a Milky Way type halo the dynamical fric- 
tion time is typically of a few Gyr, while the coalescence time for black-hole 
binaries with masses ~ 10 6 Mq (roughly the mass of the Mi lky Way MBH) 
is expected to b e < 10 7 vr jSes ana. Haardt, & Madau 20071). 

9 See however IVolonteri. Haardt, & Madau 1 2003) f° r a model that does 
not make this assumption. That model, however, only includes dark-matter 
halos and MBHs, and does not attempt to describe the formation of galaxies, 
which makes the implementation of a realistic coalescence timescale for 
MBH binaries much si mpler than in our case . 

10 We note that while Barausse & Rezzolla 12009) also present a formula 
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Because during MBH mergers a fraction of the total mass is 
radiated in gravitational waves, the mass of the MBH remnant is 
smaller than the initial mass of the binary. To account for this effect, 
we use the formula o flReisswig et alj 12009). which accurately pre- 
dicts the final mass for equal-mass black-hole binaries with spins 
aligned or anti-aligned with the orbital angular momentum: 

M bh , fln = [1 - (p + 2 Pl a + Ap 2 a)] (M bh ,i + Af bh , 2 ) (56) 

where po = 4.826x 10~ 2 ,pi = 1.559x 10~ 2 ,p 2 = 0.485xl0~ 2 
and a = (ai + ai2)/2 (ai and a 2 being the projections of the spin 
parameters on the direction L of the orbital angular momentum). 
For binaries with non -aligned spins or non-equa l masses, we use 
instead the formula of lTichv & Marron etti (2008): 

M bhjfin = (M bh ,i + M b h, 2 )x 

[1 + 4u(mo — 1) + lGmiu 2 (ai cos /3 + a 2 cos 7)] (57) 

where m = 0.9515 ± 0.001 and ra x = -0.013 ± 0.007. 

Also, because the anisotropic emission of gravitational waves 
during the merger of two generic black holes carries linear mo- 
mentum away from the system, the MBH remnant is imparted 
a recoil velocity ("kick"). There has been much controversy on 
the dependence of t his kic k velocity on the mass ratio, with 
Lousto & Zlochower d2009h advocating a scaling with v 2 , and 
Baker et alj ( 120081) finding a scaling with v 3 . It seems, how- 
ever, t hat t his discrepancy or iginated because lLousto & Zlochoweil 
(2009) and lBaker et at] (2008) considered two different regions of 
the parameter space. A compre hensive formula valid in both re- 
gions was recently proposed bv lvan Mete r etal.l ( l2010h : 

V rcco ii = v m ei + v±(cos£e! + sin£e 2 ) + v\\ e 3 , (58) 
v m = Av 2 sJ\ -4u(l + Bu), 

vj_ = H Y^- q ( a l - q<4) > 

«ij = — — l qa2 cos ( ( P 2 - * 2 ) - a i cosi^i - ®l)J 

+ [<f a 2 cos(<?!> 2 - $2) + ai cos(0i - . 

Here ei and {=3 are orthogonal unit vectors respectively in the di- 
rection of separation and along the orbital axis just before merger, 
and e 2 = ei x e.3 is a third unit vector orthogonal to them; a- 
is the projection of the spin parameter en of black hole i along 
the orbital angular momentum, while is the magnitude of its 
projection aj~ onto the orbital plane; <j>i is the angle of af" with 
respect to a reference angle representing the separation of the black 
holes, as measured at some point before the merger, while $i repre- 
sents the amount by which this angle precesses before the merger. 



predicting the orientation of the MBH remnant's spin, that information 
is not necessary in our model. As we will mention later on, in fact, our 
model assumes that the spins of the two MBHs are aligned with the an- 
gular momentum of the circumbinary disk in a gas-rich environment, in 
whic h case the merger a lways produces a final spin in the same direc- 
tion jRezzolla et al. 2008). In a gas-poor environment, instead, we assume 
that the spins of the two MBHs are randomly oriente d, in which case we 
may apply the formula of Barausse & Rezzolla 1 2009) for the final-spin di- 
rection. However, because accretion is chaotic in a gas-poor environment 
(cf. sec. |2.2.4"1 , the information about the final spin direction is not neces- 
sary in our model. 



The angles $i depend on the mass ratio and on the initial sepa- 
ration, and_jriust_be_determined with a numerical-relativity simu- 
lation ( ivan Meter et alj2010h . which seriously undermines the pre- 
dictive power of Eq. J58t . The contribution v m to the kick is dubbed 
the "mass-asymmetry contribution", because it does not depend on 
the spins and it disappears for equal mass binaries, while v± and 
Vm are the "spin contributions", which produce kicks perpendicu- 
lar and parallel to the orbital angular momentum. The angle be- 
tween the mass asymmetry and spin contributions is measured by 
£ = 215° ± 5°, while the other fitting parameters take the values 
A = 1.35 x 10 4 kms~ 1 , B = -1.48, H = 7540 ± 160 kms -1 , 
K 2 = 3.21±0.16x 10 4 kmr 1 ,^ = 1.09±0.05 x 10 5 kms^ 1 
and K s = 1.54 ±0.18 x 10 4 kms _1 . 

Because it depends on quantities measured just before the 
black-hole merger, rather than defined at large separations (as was 
e.g. the case for the formula d54t fo r the final spin, cf. the discus- 
sion in iBarausse & Rezzollal d2009h ). Eq. J58t cannot be applied 
unambiguously in our model. In fact, the orbital axis and orbital 
separation directions just before merger entering Eq. l |58l l, as well 
as the angles <£>i, can in general be determined only with a full 
numerical-relativity simulation. An exception is given by binaries 
with aligned spins, in which case the magnitude of the recoil ve- 
locity is independent of <f>i and <I>i, as can be seen from Eq. l !58t . 
and the orbital axis remains unchanged during the binary's evo- 
lution. In the genera l case, we assume (somewhat arbitrarily, cf. 
IBarausse & Rezzollal d2009h ) that the orbital axis just before the 
merger is parallel to the orbital axis at large separations, and be- 
cause 02 — 01 = X (.X being the angle between a^; and a 2 ), 
defining the phases Ai = cj>i — <E>i we can write </>i — $1 = Ai 
and 4>2 — <3? 2 = y + A 2 i n Eq. (|58[>. Using (again, a bit arb itrar- 
ily, cf. ISchnittmanl d2004l) : iKesden. Sperhake, & BertH d20ld) ) the 
value of x at large separation, we only need the phases Ai to eval- 
uate Eq. {58). In this paper, we choose these phases randomly from 
a uniform distribution! 11 ! 

The recoil velocity can be as large as 4000 km s -1 for 
an equal-mass binary with antialigned equal s pins lying on the 
equato r ial plane ("supe r kick c onfiguration", see ICampanelli et al.l 
(2007); Gonzal ez et all d2007l) ). Such a large velocity would be 
sufficient to eject the MBH remnant from the galactic nucleus. 
When the spins are aligned with one another and with the orbital 
angular momentum, instead, the recoil velocity is much smaller 
and typically insufficient to allow the remnant MBH to escape 
from the galaxy. At each black-hole merger, we therefore check 
whether the MBH remnant is retained by the spheroidal component 
of the galaxy, thereby remaining available to accrete the cold gas 
brought in by radiation drag when star formation occurs in the 
bulge [cf. Eq. l|33|ll, or whether it escapes. Therefore, after each 
black-hole merger we compare the recoil velocity Kccoii with the 
escape velocity V CS c = \/2(j>, where <j> is the potential due to the 
bulge (including both stars and gas) and the reservoir, evaluated 
at the galactic center. Because in the standard galaxy-formation 
model that we adopt here, the galactic disks form first and give 
rise to the bulges only later as a result of instabilities and major 
mergers, black-hole mergers at high redshifts will happen in 
disk-dominated galaxies, resulting in very large fractions of MBHs 



11 Another choice may be to set Ai = A 2 = 0, and we have verified 
that our results do not change significantly if we make this assumption, 
thus confirming the intuitive expectation that the astrophysical effects of 
the gravitational recoil depend more on its overall scaling with the mass 
ratio v than on the phases Ai . 
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light seeds 


heavy seeds 




0.7 


0.7 


/jet 


10 


10 


^res 


1.1 X 10" 2 


1.1 x 10~ 2 


^accr 


4.8 x 10 8 yr 


4.8 x 10 s yr 



Table 1. The calibrated values of the free parameters of the model. These 
values are used to produce the figures. 



being ejected. We will discuss this more in detail in Sec. [4] and 
hint at how this prediction may change in alternative scenarios 
of galaxy formation such a s the " tw o-phase mod e l" put forward 
in ICook. Lapi. & Granatd d2009h : ICook et alj fcOlOallbl) . As 
we will see in the next sections, however, even in our current 
"standard" model the gravitational recoil does not eject all the 
MBHs from their hosts. This is because the occupation fraction 
of black-hole seeds at high-redshifts is smaller than 1 (cf, in 
Sec. 12. II the prescriptions that we adopt to populate the halos with 
black-hole seeds at high redshifts), whic h is enough to ensure 
that M BHs survive to low redshift s ([L ippai, Fr ei. & Haimanl 



20091; Menou, Haiman, & NaravanatJ I2001I ; IVolonteril 120071 ; 



Volonteri, Gultekin, & Dottill2010h . Moreover, as emphasized by 
Schnittmanl ( 2007), even if we populated all halos with black-hole 
seeds at high redshift, the gravitational recoil in the first generation 
of mergers would automatically decrease the MBH occupation 
fraction, and even in the case of very high ejection probabilities the 
occupation fraction would settle to ~ 50 % in a few more merger 
generations. 

Also, we stress that our recipe to determine whether a MBH is 
expelled from the galaxy may be overly pessimistic. In fact, as men- 
tioned above, here we are interested in checking whether the MBH 
is retained in the spheroid, because it seems reasonable that a MBH 
wandering in the halo would neither accrete gas from the circumnu- 
clear reservoir, nor form MBH binaries after galaxy mergers. While 
this viewpoint may be sufficient for our purposes (since we want to 
study the mass and spin evolution of the nuclear MBHs and not of 
these halo MBH population), our ejection rates are clearly too large 
(especially at high z, when bulges are small) if one is interested in 
the number of MBHs that are expelled from the whole composite 
system (dark-matter halo and baryonic structures). 

In order to apply the formulas for the MBH rem- 
nant's mass, spin and kick velocity that we have just re- 
viewed, we need to specify not only the masses and spins 
of the two progenitor MBHs, but also t heir relative orienta- 
tion at large separations. It is well known teard een & Pe ttersonl 
1971 iBogdanovic, Reynolds, & Milletl l2007t iPerego et al]|2009l ; 
Dotti et alj|2010al : iMaio et alj|2012ri that if the binary's inspiral 
preceding the coalescence happens in a circumbinary disk (which 
is expected to be present in gas-rich environments), the gravito- 
magnetic torque exerted by the disk aligns the MBH spins with 
the disk's orbital angular momentum. The details of this alignment 
are weakly dependent on the equation of state of the circumbinary 
disk, with a "cold" disk (i.e. one with polytropic index T = 7/5, 
approximating a gas with solar metallicity heated by a starburst) 
resulting in a residual angle < 10° between the spins and the 
disk's angular momentum, and with a warmer disk (polytropic in- 
dex r = 5/3, corresponding to an adiabatic monatomic gas and 
thus to a scenario where radiative cooling i s suppressed duri ng the 
merger) leading to residual angles < 30° (Dotti et al. 2010jJ). As 
a result, when the dynamical interaction between the binary and 
the gas creates a low-den sity region at a separa tion of about 0.1 pc 
(the so called "gap", see I Gould & Rixl bOOOl) ). the two spins are 



aligned with one another and with the orbital angular momentum 
of the binary to within 10 — 30°. Because the gas density in the 
gap is very low, the subsequent evolution is driven by purely grav- 
itational effects which tend to further a lign the spins jSchnittmanl 
|2004 iKesden. Sperhake. & Bertill201fj|) . If the MBH merger hap- 
pens instead in a gas-poor environment, the two spins are expected 
to be randomly distributed. 

Existing models for the spin evolution of MBHs do not at- 
tempt to distinguish be tween gas-rich ("wet") an d gas-poor ("dry") 
mergers. For instance. iBerti & Volonte ri (2008) have to consider 
the two possibilities separately (i.e., either all mergers are assumed 
to be wet or all mergers are assumed to be dry) because their model 
only includes dark-matter halos and MBH s and does n ot describe 
the baryonic components; Fanidakis et al. ( 2010. [201 lh model in 
great detail the gravito-magnetic interaction of a single MBH with 
its own accretion disk, but do not include the disk's effect on the 
spin alignment prior to a black-hole merger, and consider only 
dry mergers. Because it keeps consistently track of the evolution 
of the baryonic matter, our model offers a natural way to distin- 
guish between wet and dry MBH mergers. This is important for 
the spin evolution of MBHs, because the spin distribution of the 
MBH population coming from a model with only wet mergers dif- 
fers drastically fr om that coming from a m odel with only dry merg- 
ers (see Fig. 4 in lBerti & Volonteril ( [2008l) ). Also, dry mergers tend 
to give large kick velocities, in principle sufficient to eject the rem- 
nant MBH from the galactic nucleus, because the spins prior to 
merger are randomly oriented. Wet mergers, where the spins ar e 
aligned, give much smaller recoil velocities jPotti et aljhoiOal) . 
Therefore, distinguishing between wet and dry mergers is vital to 
predict the MBH occupation number and the number of event rates 
for gravitational- wave detectors. 

More specifically, in our model we can discriminate between 
dry and wet mergers by comparing the MBH binary's mass to that 
of the circumbinary disk. This comparison needs to be performed 
carefully, however, because as explained above the MBH merger 
happens after the galactic merger, although for simplicity this de- 
lay is not implemented in our model. Major galaxy mergers are 
typically accompanied by starbursts, which in our model occur be- 
cause the galactic disks are disrupted and feed the gaseous bulges, 
where star formation is very efficient. As a result, large quantities of 
cold gas are forced into the circumbinary disk (i.e. the "reservoir" 
described in the previous sections) by radiation drag [cf. Eq. t33H , 
The circumbinary disk is therefore more massive during the MBH 
inspiral than at the time of the galactic merger. We can therefore 
approximate its mass as A/ rcs + A rcs {Mi,^ s /t s i )tg, where M les 
and M;, igas are the masses of the reservoir and bulge right after 
the galactic merger [cf. Eqs. d52b and I53H . M(, !gas /i s f is approx- 
imately the SFR in the bulge (t B f being the dynamical time of the 
gaseous bulge) and tg is the characteristic time of the MBH in- 
spiral (i.e. the free fall time of the bulge, including both stars and 
gas). During the time tg, the two MBHs continue accreting at the 
expense of the circumbinary disk, and thus at a time tg after the 
galactic merger the mass of the binary is roughly 



M binary « M bh ,i exp ( —2- ) + M bhj2 exp ( — ] (59) 



VEdd,] 



\tedd,: 



(iEdd.i and iEdd,2 being the Eddington-accretion timescales of the 
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Figur e 1. The lo cal stellar mass function (left) and baryon mass function (right) of our model compared to the observational fits of Bell et al. 1 2003^) and 
Bell etal] J2003bl) . The observational uncertainties are shown with a shaded yellow area. These results are for the light-seed scenario, but the heavy-seed one 
gives similar results. 
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Figure 2. The local gas fr action of our model compared with t he obser- 
vational parameterization of lBaldrvTG lazebrook. & Driver (2008). The ob- 
servational uncertainties are shown with a shaded yellow area. These results 
are for the light-seed scenario, but the heavy-seed one gives similar results. 



two MBHs), to be compared with a mass 

M c i r c disk ~ A/res + ^4rcs ( r^— ■ ) tff 



M b h,l 



exp 



/ tff 



V^Edd,! 



- 1 



- M b h,2 



exp 



tff 



^Edd,: 



(60) 



for the circumbinary disk. We therefore dub a merger "dry" if 
Afbinary > Af c i rc di s k, an d m tms case we assume that the spins 
of the two MBHs are randomly oriented at large separations. We 
instead consider a merger as "wet" if Mbinary < M c i rc disk, and 
in this case we assume that the spins are exactly aligned with the 
orbital angular momentum. 



3 TESTING THE MODEL AGAINST OBSERVATIONS 

Semianalytical galaxy-formation models require in general several 
free parameters to reproduce accurately a large range of observa- 



tional data (cf. for instance fBenson & Bowed ( |2010|) ). which makes 
their calibration rather involved. Our focus here, however, is not 
on the galactic properties per se, but rather on the spin evolution 
of the MBH population. We therefore allow a limited number of 
free parameters in the galaxy formation model that we described 
in the previous sections, and calibrate them against a limited but 
representative number of observations, both at z = and at higher 
redshifts. 

In particular, we allow the following four free parameters 
to vary: (i) the supernova feedback efficiency esN [cf. Eqs. d27t 
and t30H. which ranges from to 1 and which we assume to be 
the same for the feedback in the disks and in the bulges; (ii) the 
"fudge" factor /j ct [cf. Eqs. | |4U and J47H parameterizing the un- 
certainties in the Blandford-Znajeck effect (i.e. the strength of the 
magnetic field and the higher-order terms in the black-hole spin); 
this factor is assumed to vary between 0.1 and 10, and we as- 
sume it to be the same for the radio and quasar-mode feedback; 
( Hi) the normalization factor A ICS that describes the strength of the 
radiation drag regulating the growth of the circumnuclear reser- 
voir [cf. Eq. {33}]; in order to reproduce the Mbh — a relation 
and the "Magorrian" Mbh — Af& relation , it should be ^4 r ea > 
M bh /M b , and beca use M bh /M b ~ 10~ 3 (Ma gorrian et alJI 1993; 
Hari ng & Rix1l2004l) . we expect this paramete r to be on the orde r 
of 10~ 2 -10~ 3 (cf . al so iGranato et all j2004h; iLapi et al.1 d2006t) ; 
ICook et alj ( koiOallbh ; iHaiman. Ciotti. & Ostrikerl i2004l) ); (iv) the 
timescale t accr regulating the inflow of the circumnuclear reser- 
voir gas into the MBH pc-scale accretion disk [cf. Eq. {34)]; be- 
cau se the quasar bolometric luminosity p eaks at z w 2 (see 
e.g. iHopkins. Richards. & Hemquistl d2007t) ). this timescale is ex- 
pected to be on the order of the age of the universe at z ~ 2, i.e. 
taccr ~ 10 9 yr. As we will show in the next sections, these four pa- 
rameters are non-degenerate, because each of them (roughly) regu- 
lates the predictions for different observables, and this makes their 
calibration rather straightforward. 

The observables that we will consider are, at z = 0: the stel- 
lar and baryon mass functions, the MBH mass function, the gas to 
stellar mass ratio, the SFR, the Mbh — cr relation, the "Magorrian" 
Mbh — Mb relation, and the fraction of elliptical, spiral and irreg- 
ular galaxies. Also, we will test our model at z > by looking at 
the SFR density (SFRD) and at the quasar bolometric luminosity 
density as a function of z. The four free parameters of the model 
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Figure 3. Our model's SFR (at 2 = 0) compared to the observational results of Brinchmann et al. 12004). The observational 95% confidence region is 
shown with a shaded yellow area. In order to highlight possible selection effects, we show the output of our model including only the central galaxies with 
SFR> 10 -6 .Mo/yr (left panel) and both the central and satellite galaxies with SFR> lO _6 M0/yr (right panel). These results are for the light-seed 
scenario, but the heavy-seed one gives similar results. 
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Figure 4. The fraction of elliptical, spiral and irregular galaxies, at 2 = 0, 
as a function of stellar mass. The symbols are data from ConseliceJ i2006l) 
(squares: ellipticals, circles: spirals, stars: irregulars), while the predictions 
of our model are shown with thick lines (heavy seeds) and thin lines (light 
seeds). 



that we mentioned above are chosen to optimize the agreement with 
these observables, and their values are given in Tabled 

To obtain the results shown in the next sections, for each of 
our two models (light or heavy black-hole seeds) we have produced 
~ 1100 halos with masses at z = between 10 10 and 10 15 M©, 
and we have weighed the contribut ion of each one of them using the 
Sheth-Tormen halo mass function dShefh & Tormenlll999l , l2002l) at 
z = 0. We stress that each of these halos corresponds, at z = 0, 
either to a single galaxy (if there are no satellites) or to a group or 
cluster of galaxies (if satellites are present). 

3.1 Observables at z = 

For the local stellar and baryon mass func tions, we c ompare our 
results wi t h the S chechter function fits by iBell et all d2003al) and 
iBell et all J2003bl) , which were obtained with a large sample of 
galaxies from the Two Micron All Sky Survey (2MASS) and the 
Sloan Digital Sky Survey (SDSS), using simple models to convert 



the optical and near infrared galaxy luminosities into stellar masses 
and assuming a universally-applicable stellar IMF. These fits are 
shown in Fig.Q](left: stellar mass function; right: baryon mass func- 
tion) with orange thin solid lines, together with shaded yellow ar- 
eas representing the observational errors. In particular, we have as- 
sumed 0.3 dex errors in the stellar and baryon mass determinations, 
and 0.3 dex in the value of the mass functions to account for the 
statistical errors as well as the systematic uncer tainties (m a inly du e 
to the mass-to-light ratios: see discussions in IBell et all <2003al) . 



IBell et al.1 fe003bi) and also Kannapp an"& Gawisen 1 J20071) ). Also 



presented in Fig.[T]is the output of our model in the light-seed sce- 
nario (the heavy-seed scenario yields similar results), together with 
the Poissonian errors due to to finite sample of galaxies that we 
simulate. As can be seen, the model's results are generally within 
the observational errors, although the slope is significantly differ- 
ent from the observations both at the low and high-mass ends. We 
stress that the agreement of the model's mass functions with the 
observational estimates is the product of different and competing 
processes. More specifically, the ionizing radiation background and 
the supernova feedback are effective at quenching star formation 
and reducing the baryonic content in small-mass systems, while 
the quasar and radio-mode feedback, as well as the ram pressure 
and clumpy accretion, are responsible for the sharp decrease of the 
mass function at large masses, which is significantly faster than 
what would be expected from the behavior of the Sheth-Tormen 
halo mass function alone. In particular, to reproduce the high-mass 
end of the stellar and baryon mass functions it is crucial to calibrate 
the parameter /j ct , which regulates the radio and quasar mode feed- 
back in our model. 

In Fig. [2] we test our model against observational results for 
the ratio between gaseous and stellar masses ("gas fraction") at 
2 = 0. We show the comparison for the light-seed scenario, but 
the heavy-seed scenario yields simil ar results. In particular, we use 
the ga s-fraction parameterization of iBaldrv. Glazebrook. & Drivej 
(2008) (orange thin solid line), a nd we represent its observa- 
tional uncertainties (see Fig. 11 of IBaldrv. Glazebrook. & Drivej 
d2008l)) with a shaded yellow area . Because the data used by 
Baldry, Glazebr ook. & Drivej {2008) include only field galaxies, 
when calculating the gas fraction in our model we only include cen- 
tral galaxies whose stellar mass is at least 90% of the total stellar 
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Figur e 5. The local MB H mass function predicted by our model in the light-seed (left) and heavy-seed (right) scenario, compared to the observational estimate 
bv lMarconi e t al. 1 2004) The observational uncertainties are shown with a shaded yellow area. These results are for the light-seed scenario, but the heavy-seed 
one gives similar results. 



mass of the s ystem (central galaxy plus all the sate llites). Also, the 
data used by Baldrv, Glazebrook, & Driven 1 (2008) do not include 
gas-poor ellipticals, but those galaxies are only a significant field 
population at M sta , vs > 10 1 Mq, which is outside the range of 
Fig.(2] Because of the scatter of our model's predictions, we present 
its median output at a given stellar mass, as well as its 70% and 90% 
confidence-level regions, i.e. the regions containing respectively 
70% and 90% of the galaxies produced by our model at a given 
stellar mass. As can be seen, the model is in good agreement with 
the observational data. In particular, in order to obtain a good agree- 
ment with the observed gas fraction at Af s t a rs ~ 10 7 — 10 8 M©, it 
is crucial to tune the parameter esN, which regulates the efficiency 
of the supernova feedback. 

The SFR at z = as a function of stellar mass, as ob- 
tained fro m a survey ~ 10 5 gala xies with measurable SFR in 
the SDSS (Brinchma nn et alj|2004) . is shown in Fig. [3] The or- 
ange thin solid line represents the median SFR observed at a given 
stellar mass, while the shaded yellow area represents the obser- 
vational 95% confidence level (i.e. the region containing 95% of 
the galaxies at a given stellar mass.). The median SFR predicted 
by our model at a given stellar mass (in the light-seed scenario: 
the heavy-seed one gives again similar results) is represented by 
the thick solid blue line, while the thin red dashed and thin purple 
dot-dashed lines represent the model's 70% and 90% confidence 
re gions. In particular , becau se the median and confidence region 
of Brinchmann et 

Id] 12004) may depend on the composition of 
their sample of galaxies, we have plotted the output of our model 
when considering only the central galaxies, where star formation 
is more intense (left panel), and when considering both the central 
galaxies and the satellites (right panel). In both cases, we have ne- 
gl ected galaxies with SFR le ss than lO~ 6 A/0 /yr, since the sample 
of iBrinchmann e t al. (2004) only includes galaxies with measur- 
able SFR. 

As can be seen, irrespective of which of the two samples we 
consider, the median predictions of our model lie within the 95%- 
confidence region of the observational data, and reproduce the flat- 
tening of the SFR at high masses (due to the combined effect of 
quasar and radio-mode feedback, clumpy accretion and ram pres- 
sure). However, our predictions are significantly lower than the me- 
dian of the observations, especially at small stellar masses, essen- 
tially because they present a steeper slope. More specifically, our 



model predicts SFR oc M s t ars , while the data by Brinchman n et al.l 
d2004|) suggest SFR oc M s " ars , with n w 0.7. While other ob- 
servati onal data point at a slope n slightly l ess than 1 (for in- 
stance. lElbaz et all J2007h find n s a 0.77 and lSalim etal] J2007h 
find n sa 0.65). lElbaz et all < l201ll) find that at z = the SFR- 
M s tars relation is well fitted by SFR= Af stars /(4 x 10 9 yr) for 
9.5 < log 10 A/ S tars < 11.5. This parameterization is represented 
in Fig. [3] with a dashed thick black line. Also, there is some evi- 
dence t hat the slope n is ra ther sensitive to selection effects: for in- 
stance, iKarim et al.l d201ll) find that for highly active star-forming 
galaxies the exponent n approaches 1, although it is not clear 
whether these galaxies are representative of the entire star-forming 
population. In general, the slope n ~ 1 of our model's predictions 
is due to our star formation prescriptions of Sec. 12. 2. 3"! and is there- 
fore commo n to most semianalytical galaxy-formatio n models (see 
for instance iDutton. van den Bosch. & Dekell feoid) ). While it is 
in principle possible t o change our st ar formation prescription to 
obtain a milder slope JShi et al]|201lh . at this point it is not clear 
whether this is needed, because it seems that the issue is still not 
settled from the observational point of view. 

In Fig. [4] we plot the fraction of elliptical, spiral and irregular 
gala xies as a function of stellar mass. The symbols are actual data 
from lConsellcel (12006) (squares: ellipticals, circles: spirals, stars: ir- 
regulars), while the predictions of our model are shown with thick 
lines (heavy seeds) and thin l ines (light seeds) . More specifically, 
within our model we follow iGuo et alJ d201ll) and use the ratio 
Mb/Mtot (with Aftot = Mb + Md) to discriminate different mor- 
phologies. In particular, we identify galaxies with Mb/Mtot > 0.7 
with ellipticals and represent them with solid lines; galaxies with 
0.03 < Mb/Mtot < 0.7 with spirals and represent them with 
dashed lines; galaxies with Mb/Mtot < 0.03 with ext reme late- 
type o r pure-disk galaxies (and therefore with irregulars dGuo et all 
201 1 )). and represent them with dotted lines. Remarkably, in spite 
of this simplistic classification, our model seems to reproduce the 
observed morphological fractions, at least qualitatively. 

To test our predictions for the MBH population at z = 0, 
we look at the local MBH mass function, at the Afbh — a rela- 
tion between the black-h ole mass and the line-of-sight velocity dis- 
persion a of the bul ge (Ferrarese & Merritt 2000; Gebhard t et al. 
200d:lGultekin et al.ll2009l ; lGraharrJl2008l : lGraham et alj|201 lhlHvj 



2008), and at the relation between the MBH mass and the bulge 
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Figur e 6. The predic tion s of our model, in the light-seed scenario, for the Mbh — a and M^h — Mj, relations at z = 0, compared to the observational fits of 
Giiltekin et al. (2009) and Haring & Rix (2004), respectively. The observational uncertainties are shown with a shaded yellow area. 



9.5 

9 
8.5 

8 
7.5 

7 
6.5 

6 



Best fit to data 
Model's median 
70% confidence 
90% confidence 




2.1 2.2 2.3 
log 10 a [ km/s ] 



2.5 




j t. i_ 



9.5 10 10.5 11 11.5 12 12.5 

lo 9l0 M bulge I M J 



Figure 7. The same as in Fig. [6] but for the heavy-seed scenario. 



mass, ini tially proposed by Magorria n et al.l dl998h and later up- 
dated by Hari ng & Rixl d2004l) . The predictions of our model for 
these observables mainly depend on the normalization factor A Tes 
regulating the growth of the circumnuclear reservoir and therefore 
of the MBHs, and we calibrate this parameter to obtain agree- 
ment with t he observation s . In p articular, in Fig. [5] we show the 
estimate by iMarconi et al. I d2004l) for the MBH mass function at 
z — with a thin solid orange line, and we assume observa- 
tional uncertainties, represented by shaded yellow areas, of 0.3 
dex for the MBH mass and 0.4 dex for th e mass function (cf. 
IShankar. Weinberg. & Miralda- Escudd d2009h . Fig. 5). The predic- 
tions of our model are instead shown in the left panel for the 
light-seed scenario, and in the right-panel for the heavy-seed sce- 
nario, with error bars representing the Poissonian errors due to to 
finite sample of galaxies that we simulate. As can be seen, both 
scenarios produce MBH mass functions that are compatible with 
the observational estimate, but they yield different predictions for 
M bh < 1O 6 M , where the mass function is unconstrained by ob- 
servations. In particular, the heavy-seed scenario obviously predicts 
a large number of MBH with M b h > M soed = 10 5 M Q , while the 
light-seed scenario predicts a flatter mass function. 

In the left panel of Figs. [6] (light-seed scenario) and[7] (heavy- 
seed scenario), we show the parameterization of iGultekin et al.l 



1 2009) for the Mbh — o relation orig i nally discovered by 



iFerrarese & Merritj d2000l) : iGebhardt et all d2000h, while in the 
right panel we show the parameterization of Haring & Rix ( 200' 
for the Mbh — Mj, relation, originally suggested bv lMagorrian et 
dl998h . In particular, both parameterizations are represented by thin 
solid orange lines, while the lcr observational scatter is denoted by 
a yellow shaded area. We also show the median of our model's pre- 
dictions for Mbh at any given a (Mb), as well as their 70% and 
90% confidence regions. 

More specifically, in order to calculate the median and confi- 
dence regions in our model we have neglected the MBHs residing 
in satellite galaxies. Also, even for the central galaxies, we have 
only considered the MBHs residing in elliptical galax ies (which 
we id entify again with ones having Mb/M to t > 0.7 dGuo et all 
l201lh ). This is because the Mbh — o relation is kno wn to present a 



large s catter for late-type galaxies (cf. discu s sion in Giiltekin et al 
(2009) and their Fig. 1; see also iGraharr] J2008h : iGraham et al. 
( 201 lh : [HuI J2008h'). and the M b h - M b relation has been estab- 
lished in lHaring & Rixl (2004) using a sample of 30 mostly ellip- 
tical galaxies. The predictions of our model are obtained using the 
virial theorem to relate the bulge velocity dispersion V^ u i ge to the 
bulge gravitational potential (which we calculate from the bulge 
density d21t). and then calculating the line-of-site velocity disper- 



© 0000 RAS, MNRAS 000, 000-000 



The evolution of massive black holes and their spins in their galactic hosts 17 



sion as a 2 = V^ ulge /K, where the correction factor K depends 
on the type of orbits of the stars in the bulge. For instance, if the 
orbits were mainly along the line of sight, one would have K w 1, 
but K ~ 3 for nearly circular and isotropic orbits. Here, follow- 
ma iRood Jl970h , we assume K — 2.1. Overall, as can be seen, 
both the light-seed and heavy-seed scenario reproduce the obser- 
vational relations only qualitatively, because they predict a signif- 
icant number of "outliers" below the observed correlations, with 
masses Mbh < 10 8 . These outliers become even more numerous 
if one includes the spiral and irregulars galaxies in the analysis, or 
if one also accounts for the MBHs in satellite galaxies. Physically, 
they therefore represent MBHs that will settle on the Mbh — a and 
Mbh — Mb relations in the future, after merging with a MBH re- 
siding in a central galaxy already on the Mbh — cr and Mbh — Mb 
relations, or after a burst of star formation in the bulge following a 
disk disruption or a major merger. This pop ulation of outliers is als o 
found for instance in the model of Volont eri & Nataraianl J2009t) . 
but remarka bly not in the "two-ph a se" g a laxy-formati on model 



proposed bv lCook. Lapi, & Granatd j2009l):|Cook et all fcoiodlbl) 
(see in particular Fig. 7 of lCook et alj yOlObJ)), where bulges form 
at high redshifts and MBHs settle on the M b h — cr and M hh — M b 
relations earlier. We will discuss this "two-phase" model and its 
possible consequences on the MBH redshift evolution more in de- 
tail in Sec.|4] 



3.2 Observables at z > 

In this section we will examine the output of our model for two 
observables at z > 0, namely: the SFR density (SFRD) and the 
bolometric quasar luminosity density. 

There are at least two classes of methods to determine the 
redshift evolution of the SFRD, also known as the st ar-formation 
"cosmic history" ( Mada u et alJl996l : lLillv et all 1996). On the one 
hand, the SFRD can be measured directly at various redshifts using 
instantaneous indicators. More specifically, emission from massive 
stars, which are short-lived compared to the typical star-formation 
time scales, can be u s ed to e xtrapolate the IMF and obta i n the SFR 
(see iHopkinj d2004 l2007h : iHopkins & Beacoml < |2006|. l2008h for 
compilations of recent data). A similar procedure ca n be followed 
with c ore-collapse supernova rates (see for instance iDahlen et al.l 
HH)), for these stars too are short-lived compared to star forma- 
tion. Alternatively, one can measure the distribution of stellar ages 
in nearby galaxies and use stellar population synt hesis models to re - 
construct the star-formation cosmic history (see lPanter et alj fc007l) 
for a recent analysis). 

On the other hand, integrating the cosmic star-formation his- 
tory over redshift and correcting for the mass loss through super- 
novae and stellar winds yields the stellar-mass density as a function 
of redshift. Conversely, the cosmic star-formation history can be ex- 
tracted from the stellar-mass density evolution, which can be mea- 
sured independently with galaxy surveys and is sensitive to a larger 
range of masses than instantaneous indicators, which only probe 
the most massive stars. Moreover, instantaneous indicators are sub- 
ject to greater uncertainties due to the effects of dust obscuration. A 
compilation of recent measurement s of the stellar-mass density as a 
functio n of redshift is presented bv lWilkins. Trentham. & Hopkins! 
(2008), who showed that the inferred star-formation history agrees 
with the one derived from instantaneous indic ators for z < 0.7, at 
least f or suitable choices of the IMF (see also Hopki ns"& Beacoml 
(2008)). At higher redshift, however, the instantaneous indicators 
give larger SFRDs than the evolution of the stellar-mass density, 
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Figure 8. The cosmic st ar-formation history predicted by our model, com- 
pared to that derived by Wilkin s. Trentham, & Ho pkins ( 2008) from obser- 
vations of the stellar mass-density and assuming an evolving IMF (yellow 
and orange shaded areas, representing the lcr and 3<r uncertainty regions 
of t he observations). Also shown are fits to instantaneous SFRD indica- 
tors iwilkins. TrenthamT & Hopkins 2008), assuming either a universal or 
an evolving IMF. 
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Figure 9. The prediction of our model for the bolometric quasar 
luminosity density as a function o f redshift, compared to data from 
Hopkins, Richards, & Hernquist |2003). 



and thi s discrepancy peaks at z m 3, where the difference is ~ 0.6 
dex iWilkins. Trentham. & Hopkinsll2008l : lwilkins et all2008l) . 

There are various possible explanations for this discrepancy, 
namely uncertainties in the effect of dust on the stel lar-mass and 
SFR estimates at high-redshifts toriver et alj |2007|). or incom- 
pleteness in the measured stellar-mass density i Reddv & Steidell 
2009). Another possibility is that the IMF might evolve with 
redshift JWilkins. Trentham. & Hopkinsll2008l ; IWilkins et all2008t 



lDavejr2008l) , and there are in fact th eoretical arguments and in 
direct observational evid ence for that ivan Dokkum| [2008 : Larson 
l2005l;IWang & Daill201l1) . In Fig.[U the yellow and orange shaded 
areas represent the la and 3<r un certainty regions of the cosmic 
star-fo rmation history derived by IWilkins. Trentham. & Hopkins! 
d2008l) from the evolution of the stellar mass-density, as- 
suming an evolving IMF. The d ashed and solid lines are 
IWilkins. Trentham. & Hopkins! |2008)'s fits to instantaneous SFRD 
indicators, assuming respectively a universal and an evolving IMF. 
In this figure, we also show the output of our model. As can be 
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seen the light-seed and heavy-seed scenarios predict essentially the 
same SFRD, but these predictions agree only qualitatively with the 
observations. While it should be possible to amend our model to 
reproduce the observational data more closely, that would proba- 
bly require including more free parameters than the four that we 
consider in this paper (as of now, we have no parameters regulat- 
ing the cosmic star-formation history directly, besides those that we 
fixed by comparing to the z = observations of the previous sec- 
tion). We deem such a refinement premature due the discrepancies 
between different indicators and the problems in interpreting the 
SFRD observations that we mentioned above. 

Finally, in Fig. [9] we show data from 
Hopki ns! Richards. & Hernquistl J2007b for the bolometric 
quasar luminosity density as a function of redshift, and the 
corresponding output of our model. As can be seen both the light- 
seed and heavy-seed scenarios do a good job at reproducing the 
observations, and in particular the peak of the quasar luminosity 
at 2 w 2. We stress that while the normalization of the quasar 
luminosity density depends on the parameter ^rcs, which also 
regulates the normalization of the MBH mass function and that 
of the Afbh — a and Mbh — Mt, relations (see previous section), 
its shape crucially depends on our free parameter / aC cr, which we 
calibrate to reproduce the observations. Also, the predictions of 
our model as shown in Fig.|9]only consider MBHs with bolometric 
luminosities L > 1O 1O L . This is because lower luminosities 
cannot be observationally resolved as confusion with normal 
star-forming and starburst galaxie s becomes an issue, and therefore 
they d o not enter in the analysis o fHopki ns, Richards, & Hernquistl 
d2007h . 



4 THE CHARACTER OF MBH MERGERS 

In this section we examine the predictions of our model for the 
MBH mergers, focusing in particular on their character, i.e. whether 
they happen in gas-rich ("wet") nuclear environments (where the 
gravito-magnetic torques align the MBH spins prior to the merger) 
or in gas-poor ("dry") ones, where the spins prior to the merger are 
randomly oriented. 

In Figs.[lO](light-seed scenario) andQTJ(heavy-seed scenario) 
we present predictions for the number of mergers observed in 1 yr 
at z = 0, for the fraction of MBH remnants that are ejected from 
galactic spheroids as a result of the gravitational recoil, and for 
the fraction of wet mergers, in different ranges of the MBH binary's 
mass Afbin = Mbh i + Mbh 2 and as a function of redshift. The 
predictions for the number of mergers, however, might be regarded 
as lower limits to the actual rates, because of the prescriptions de- 
scribed in Sec. 12.11 and aimed at keeping the computational time 
needed to follow the merger trees up to high redshifts to an accept- 
able level. 

In spite of this note of caution, the results for the light-seed 
scenario seem to confirm that LISA/SGO or a similar European 
mission such as eLISA/NGO should detect at least a few MBH 
merger events during its lifetime, although a detailed analysis will 
be necessary when the details of the mission (e.g. its duration and 
sensitivity curve) are finalized. In fact, focusing for the moment on 
the 1O 4 M < M bin < 1O 6 M and M bin > 1O 6 M mass ranges 

12 As mentioned in Sec. |2. 3. II our results for the ejection rates are overly 
pessimistic if one is instead interested in the fraction of MBHs that are 
ejected from the whole composite system (dark-matter halo and baryonic 
structures). 



(which are the most relevant for LISA/SGO and eLISA/NGO), we 
notice that the merger rate s shown in F ig.llOlare qualitatively sim- 
ilar to the predictions of ISesanaT Volon teri. & Haardt {200% for 
the "BVRhf" model (cf. their Fig. 1), which yields event rates 
of a few per year for the original LISA mission. In the heavy- 
seed scenario (Fig, II U. the events in the mass range 1O 4 M < 
Mbin < 1O 6 M are instead much more numerous (~ 200 per 
year for z < 10). This is ind eed in agreement with the analysis of 
ISesana. Volonteri. & Haardl d2007l) . who found that models with 
heavy seeds and large initial seed occupation numbers (i.e. their 
"KBD" model) should give event rates of hundreds per year for 
LISA/SGO or eLISA/NGO. As for the low-mass range M bln < 
1O 4 M , our light-seed scenario predicts about 25 mergers per year 
for z < 10. (Clearly, mergers in this mass range are not present 
in the heavy-seed scenario.) These mergers may give a significant 
event rate for future third generation gravitational-wave detectors 
in the 0.1-1 Hz frequency band, such as DECIGO or t he Einstein 
Telescope dSesana et all 120091 ; iGair et al.l [20091 1201 ll) . or might 
even be marginally detectable with LISA/SGO or eLISA/NGO at 
high redshifts if they have M hin < 1O 4 M . 

As for the fraction of MBHs ejected from galactic bulges and 
for the fraction of wet mergers, Figs. UOIandll llconfirm one's in- 
tuitive expectations. In the light-seed scenario, the mergers in the 
Mbin < 1O 4 M mass range almost always result in the MBH 
remnant being ejected. This is because the mass of MBHs remains 
small (i.e. ~ M see d) only if the host galaxies are disk-dominated, 
in which case little or no star formation happens in the bulge, and 
no cold gas becomes available for the MBHs to grow. Naturally, if 
the bulge components are small they can hardly retain the MBHs. 
Also, because the MBH seeds all have the same mass, mergers be- 
tween seeds have mass ratio q ~ 1, which gives larger recoil ve- 
locities (cf. the mass-ratio dependence of Eq. J58» . Mergers in the 
1O 4 M < Mbin < 1O 6 M and M bin > 1O 6 M mass ranges, 
instead, are less likely to result in MBH ejections in the light- 
seed scenario. This is because if the MBHs have managed to grow 
beyond mere seeds, they have done so by accreting the cold gas 
brought to galactic nuclei by star formation in the bulges (via ra- 
diation drag). As a result, the bulges are more massive than in the 
case of mergers of MBH seeds, and they are more likely to retain 
the MBH remnants resulting from mergers. Also, because the MBH 
have grown to masses far from that of their seeds, the mergers are 
likely to involve mass ratios significantly different from q = 1, This 
is indeed shown in Fig.[T2] where we show the different mass ratios 
occurring in MBH mergers in the light-seed scenario, in different 
mass ranges. As can be seen, for Mbi n < 1O 4 M the coalesc- 
ing black holes often have comparable masses (essentially because 
many of them are still seed black holes, or have grown little away 
from them), while the mass ratios become more varied at higher 
masses. 

Similarly, as can be seen from Fig.QI)] the mergers in the mass 
range Afbin < lO 4 Af are wet only in ~ 30 — 40% of cases in 
the light-seed scenario, because if a significant amount of gas were 
present in the galactic center, the black holes would have rapidly 
grown beyond 1O 4 M . In the 1O 4 M < M bln < lO 6 A/ and 
Afbin > 10 6 Mq mass ranges, the mergers are mainly wet at high 
redshifts, where a lot of cold gas is present in galactic nuclei, but 
the fraction of wet mergers decreases as the cosmic evolution pro- 
gresses, because the amount of gas shrinks as a result of both ac- 
cretion by the MBH and its quasar and radio-mode feedback on the 
galaxy (as well as a result of ram pressure and clumpy accretion, at 
large halo masses and low redshifts). Such a decrease in the fraction 
of gas-rich MBH mergers with redshift is a prediction that can be in 
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Figure 10. The predictions of our model (in the light-seed scenario) for the number of MBH mergers per year of observation, for the fraction of mergers 
producing a MBH that is ejected from the galaxy as a result of the gravitational recoil, and for the fraction of gas-rich ("wet") MBH mergers, as a function of 
redshift and in different mass ranges. We notice that the number of mergers with Afbin > 10 6 Mq drops to zero for z > 10, hence the fraction of wet mergers 
and that of ejected MBHs are not defined for z > 10 in that mass range. 
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Figure 11. The same as in Fig. 1 101 but in the heavy-seed scenario. We notice that no mergers with Afbin < 10 4 A/q are present, because of the seed model. 



principle tested with gravitational-wave missions (e.g. LISA/SGO 
or eLISA/NGO). As already mentioned, the MBHs spins tend to 
be at least partially aligned in gas-rich environments due to the 
gravito-magnetic torques exerted by the gas, while the spins are 
expected to be randomly oriented and precess more strongly in in 
gas-poor ones. Missions like LISA/SGO or eLISA/NGO should be 
able to tell a binary with nearly aligned spins from a binary with 
strongly precessing ones by loo king at the higher-order harmonics 
of the gravitational waveforms jLang & Hughedl200dl200ll2008l ; 
lLang. Hughes. & Cornishll201ll) . 

The situation is similar in the heavy-seed scenario, as can be 
seen from Figs.[TT]and[T3] The mergers in the lO 4 A-/0 < Mbi n < 
1O 6 M0 mass range happen between MBHs that have barely grown 
away from their original seeds of mass M aeed = 10 5 M©. There- 
fore, these mergers usually present comparable masses, happen in 
dry environments, and likely eject the resulting MBH from the 
galactic spheroid. The mergers in the A/bi n > 1O 6 M0 mass range, 
instead, consist of MBHs that have grown significantly larger than 
their seeds thanks to a gas rich environment and to significant 
star formation in the bulge, and therefore tend to be retained in 
spheroids after mergers, and to merge in "wet" environments (at 
least at high redshifts, where circumbinary disks have not yet been 
destroyed by accretion, AGN feedback, ram pressure and clumpy 
accretion). 

We stress that in spite of the rather high ejection probability 
in the mergers between seeds that one can observe in Figs. [10] 
and [TT] our model correctly reproduces the property of MBHs 
at low redshifts, as we showed in the previous section. In other 
words, the gravitational recoil does not succeed at rooting out all 
the MBHs from their host galaxies. This is essentially because 



the occupation fraction of black-hole seeds at high-redshifts is 
smaller than 1 (cf., in Sec. 12.11 how we populate high-redshift 
halos with black-hole seeds), and this is enough to ensure 
that M BHs survive to low redshift s |Lippai, Frei, & Haiman 
20091; iMenou. Haiman. & Narayanan! l200ll ; IVolonteril 120071: 
Volon teri. Gultekin. & Dottill2010i) . Also, as noted bv lSchnittmanl 
( 2007), even if we populated all halos with black-hole seeds at high 
redshift, the gravitational recoil in the first generation of mergers 
would automatically decrease the MBH occupation fraction, and 
even in the case of very high ejection probabilities the occupation 
fraction would settle to ~ 50 % in a few more merger generations. 

Ano ther possibility would be to a d opt the galaxy-form ation 
model of ICook. Lapi. & Granatol fc009l) : ICook et alj j2010al lbh. in 
which the baryonic evolution is driven by the two-phase stru c- 
tural evolution of the dark-matter haloes dZhao et all l2003al fbl 
iMo & Maoll2004lDiemand. Kuhlen. & Madaull2007l) . and presents 
two distinct phases: an early "fast collapse" phase, where the dark- 
matter core structure is built through a series of violent merger 
events, corresponding to an epoch where baryonic material col- 
lapses to directly form spheroidal structures (bulges); and a late 
"slow collapse" phase, where potentially large amounts of ma- 
terial are added to the halo outskirts without affecting the cen- 
tral regions, giving rise to the quiescent growth of disk struc- 
tures around the previously-formed bulges. This model can po- 
tentially reproduce the observed downsizing of baryonic struc- 
tures mo re naturally than standard sem ianalytical galaxy-formation 
models dCook. Lapi. & Granatol l2009h . and may potentially have 
effects on the predictions for the MBH merger rates and mass 
evolution. As already mentioned in Sec. [3] the standard galaxy- 
formation model that we adopt here, in which disk galaxies form 
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Figure 12. The predictions of our model (in the light-seed scenario) for the distribution of mass ratios q = Mbh s/^hh 1 (where Mbh,2 < ^bh,l) m MBH 
mergers, as a function of redshift and in different mass ranges. 
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Figure 13. The same as in Fig. 1 121 but in the heavy-seed scenario. 



first and give rise to spheroids by instabilities and major merg- 
ers, predicts the presence of a significant number of outliers in 
the Afbh — a and Mbh — Mb relations, while MBHs settle on the 
Mbh — & and Mb h — Mb relations ea rlier in the case of this "two- 
phase" model (see lCook et al.l l l2010bl )). This is because spheroidal 
structures form first, and as a result of radiation drag they feed the 
MBHs at higher redshifts than in the "standard" model. Clearly, in 
such a scenario one does not have high-redshift mergers between 
galaxies with little or no bulges, and therefore the ejection rate of 
black-hole seeds should be significantly reduced. Also, the earlier 
growth of the MBHs may boost the event rates for LISA/SGO or 
eLISA/NGO in the light-seed scenario. We will explore in detail 
the effects of this alternative galaxy-formation model on the MBH 
mergers and evolution in a future paper. 

Finally, as mentioned in the introduction, the results of this 
section depend on our assumption that the MBH spins gets aligned 
by gravitomagnetic torques in gas-rich environments. However, due 
to processe s such as star formation and feedbac k in the circumnu- 
clear disk footti et alj|2010al ; iMaio et alj|2012h or clump forma- 
tion i n high-2 disk galaxies dBournaud et alj 1201 1\ ; iDubois et al.l 
l201lh . this alignment is likely to be only partial, which may result 
in higher ejection rates for the MBHs resulting from mergers. 



5 THE EVOLUTION OF THE MBH SPINS 

As a final application of our model, we study the redshift 
evolution of the MBH spins. These predictions will be read- 
ily testable by LISA/SGO or eLISA/NGO, which will mea- 
sure the black-hole masses and spins with astonishing ac- 



cur acy 



10" 



for the masses and 10 for the 



spins, 



;|Berti, Buonanno. & Willi d2005h:lLang & Hughel {2006, 2007, 
l2008h : lLang. Hughes. & Cornishl fcOllh T and without the system- 



atic uncertainties typically affecting electromagnetic (e.g. X-ray) 
determinations. 

In Figs.[l4](light-seed scenario) and[T5](heavy-seed scenario), 
we present results for the distribution of masses and spins of the 
MBHs residing in isolated galaxies or in the central galaxies of 
groups or clusters (i.e., we do not consider the MBHs residing 
in satellite galaxies), at redshifts ranging from z = 7 to z — 0. 
The color code represents the log 10 of the density of MBHs per 
unit (logarithmic) mass and unit spin, log lo (d0bh[Mpc _,:! ]/da) = 
log 10 (d 2 nbh[Mpc _3 ]/(d log 10 Mbh[MQ] da))) . As can be seen, 
in the light-seed scenario, already at z — 7 the MBH distribution 
has been skewed towards large spins from the initial uniform spin 
distribution of the seeds (still visible at M scac \ = 150Mq). This 
is because at high redshifts, where the AGN feedback is still inef- 
fective and the MBHs small, large amounts of gas are present in 
galactic nuclei, and the MBH spins grow as a result of wet, spin- 
aligned mergers (cf. Sec. [4]l and most importantly because accre- 
tion onto the MBHs is coherent and spins them up (cf. Sec. |2.2.4] l. 
We stress, as already mentioned in the introduction, that effects 
such as star formation and feedback in the circumnuclear disk, 
or the formation of clumps in high-redshift disk galaxies fed by 
cold streams, may at least partially randomize the accretion flows 
in gas-rich environments. As a result, the large spins eibh ~ 1 
shown in Figs. [14] and [15] at high redshift may be substantially re- 
duced (e.g. Mai o et al. 1 120121) show that for quasars in merger rem- 
nants, sustained accretion results asymptotically in spins parame- 
ters a bh ~ 0.7 - 0.9). 

At smaller redshifts this trends gets modified because the cold 
gas in the nuclear regions of galaxies becomes scarcer, hence merg- 
ers tend to happen in dry environments (cf. Sec. [4]l and accre- 
tion turns chaotic. Chaotic accretion, in particular, appears to be 
the main driving force behind the spin evolution in this phase, as 
can be seen from the appearance of a large number of MBHs with 
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Figure 14. The evolution of the MBH masses and spins with redshift, as predicted by our model in the light-seed scenario. The color code represents the log 10 
of the density of MBHs per unit (logarithmic) mass and unit spin, i.e. log 1{ )(d</>bh [Mpc - 3 ]/da) = log 1() (d 2 ribh[Mpc — 3 ]/(d log 10 M^[Mq] da))) . 



spin parameter abh ~ 0.1. This is indeed what would be expected 
in a chaotic-accretion scenario, wher e the black-hole spi n oscil- 
lates around a small non-zero value ( iKing & Pringiel l2006). In our 
model, the value abh ~ 0.1 is easily explained. As mentioned in 
Sec. 12.2.41 accretion turns chaotic when the mass of the gaseous 
reservoir drops below the black hole's mass ("dry" environment). 
Assuming that the MBH is almost maximally spinning as a result of 



the previous phase of coherent accretion, we can calculate the spin 
of the MBH when the reservoir has been completed accreted by in- 
tegrating Eq. l |39l > from an initial spin abh = 1, and assuming that 
the MBH accretes a mass of gas M TCB = M b h- Doing so, one gets 
a final spin abh ~ 0.14, which explains the large number of MBH 
with spin abh ~ 0.1 at low redshifts. The evolution of the spins is 
qualitatively similar in the heavy-seed scenario, with the difference 
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Figure 15. The same as in Fig. 1 141 but in the heavy-seed scenario. 



that no MBHs with Mhh < M BCC d = 1O 5 M0 are present, because 
of the seed-model described in Sec. 12. II 



We stress that our results, and in particular the dichotomy 
between almost maximal spins at high redshifts and small spins 
at z ~ 0, are qualitatively independent of our assumption that 
the seeds are initially assigned spin parameters drawn from a uni- 
form distribution —1 < abh < 1, at least for MBH masses 
A^bh ^ 3M sec d . While it is unclear whether such a spin-parameter 



distribution makes sense physically, because little is known about 
the spins of the seeds, at high redshifts the MBHs accrete coher- 
ently, and they lose memory of their initial spin after accreting a 
mass comparable to their own (i.e., if a black hole of mass Afbh ac- 
cretes coh erently, its sp in becomes maximal after accreting a mass 
< 2M bh dBardeenlll97(j|) ). 

In both the light and heavy-seed scenario, it would seem 
that the paucity of MBHs with large spins might be in contrast 
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with the iron Ka measurements of the MBH spins in MCG- 
6-30-15 and in NGC 3783, which were claimed t o be respec- 
tivel y Qbh > 0.987 (Brenneman & Reynolds 2006) and abh > 
0.9 jRevnolds et~afll201ll ; iBrenneman et al.l 1201 lh at 90% confi- 
dence level. It should be noted, howe ver, that both these mea- 
surements are still controversial, with IPatrick et alj J201ll) find- 
ing a bh = 0-49±S;?a for MCG-6-30-15 and o bh < -0.04 for 
NGC3783, and with other iron Ka measurements of MBH spins 
giving smaller values, e.g. a spin between 0.3 and 0.77, according to 
the m easurement, for Fairall 9 ( Patric ket alj|201 it iRevnolds et al.l 
2011). Even more importantly, these spin measurements are nec- 
essarily biased by selection effects, because large spins correspond 
to higher em ission efficiencies ?7(abh) and therefore high er AGN 
luminosities (Reynol ds et alj|201lt IBrenneman et alj|201ll) . More- 
over, iron Ka measurements are of course only possible in systems 
with accretion disks in the first place (i.e in AGNs), and those sys- 
tems are expected to host MBHs with high spins because coherent 
accretion spins them up to the maximal limit (cf. Sec. 12.2.41 . 



6 CONCLUSIONS 

We have utilized a semianalytical galaxy-formation model to track 
the evolution and mergers of the dark-matter halos, the IGM, the 
baryonic structures (galactic disks and spheroids, in both their 
gaseous and stellar components), and the MBHs that are thought 
to reside in the center of galaxies. The evolution of the MBHs is 
deeply entangled with that of their host galaxies, because it is the 
star formation in the galactic spheroid that funnels gas to the galac- 
tic center via e.g. radiation drag. This creates a reservoir that feeds 
the MBH, but the character of the accretion process (coherent or 
chaotic) depends on the amount of gas present in this reservoir. 
Similarly, when two MBHs merge after their host galaxies have co- 
alesced, their spins get aligned prior to the merger due to gravito- 
magnetic torques if enough gas is present in the galactic nucleus, 
while the orientation of the spins remains essentially isotropic in 
gas-poor environments. A further complication to this picture is 
that the MBHs also backreact on the larger-scale galactic evolu- 
tion, through the so-called AGN feedback, i.e. they are thought to 
quench star formation in high-mass systems by injecting energy 
into the IGM via strong jets or accretion-disk winds. Indeed, the 
AGN feedback is a crucial ingredient of modern galaxy-formation 
models, and is needed to explain the "anti-hierarchical" evolution 
(or "downsizing") of baryonic structures, i.e. the fact that the most 
massive galaxies are dominated by old stellar populations, while 
low-mass galaxies generally present young stellar populations and 
longer-lasting star-formation activity. 

In this paper, we have made an attempt to study how this 
complicated interdependence between MBHs and their host galax- 
ies affects the black-hole mass and spin evolution, considering 
both a scenario where MBHs form from "light" 150Mq seeds at 
2 ~ 15 — 20 and one where they form from "heavy" ~ 1O J M0 
seeds at z ~ 10 — 15. Besides confirming that these two scenarios 
may be observationally distinguishable with LISA/SGO or a simi- 
lar European gravitational-wave mission (eLISA/NGO) by simply 
looking at the observed event rate for MBH mergers, we have stud- 
ied the MBH mass and spin evolution in detail. In particular, we 
have determined that accretion is mostly coherent and MBH merg- 
ers happen in gas-rich environments at high redshifts, while at low 
redshifts, when AGN feedback, ram pressure and clumpy accretion 
have "sterilized" the galaxy, accretion becomes mainly chaotic and 
mergers happen in gas-poor environments. This results in a spin 



distribution that is skewed towards large spins at high redshifts, and 
towards spins abh ~ 0.1 at low redshifts, a prediction that will be 
readily testable by LISA or a similar mission. 

In principle, LISA/SGO or eLISA/NGO will also be capa- 
ble of testing our predictions for the character of MBH mergers 
as function of redshift directly, because gas-rich mergers tend to 
present aligned spins, while gas-poor mergers tend to present ran- 
domly oriented spins, and these different orientations will pro- 
duce an effect on the higher-order harmonics of the gravita- 
tional waveforms. We will examine this more in detail in future 
work aiming at calculating the eLISA/NGO signal-to-noise ra- 
tios and the effect of higher-order harmonics. In particular, to ac- 
curately model the gravitational waveforms for spinning black- 
hole binaries we will employ the effective-one-body model of 
iBarausse & Buonannol J2010l 1201 ll) . which successfully repro- 
duces the exact (i.e. nu merical-relativity) waveforms both in the 



extreme-mass ratio limit jYunes etaf 201 it Barausse et Hi] 1 20 111) 
and in the comparable-mass case jPanetall201lh . 

Finally we have briefly hinted at how the mass evolution of 
MBHs and the event rates for LISA/SGO or eLISA/NGO may 
be different in alterna tive galaxy-formation models such as the 
"two-ph ase" model of ICook. Lapi. & Granatol d2009h : ICook et all 
(2010a b). This model reproduces the downsizing of cosmic struc- 
tures more naturally than the standard galaxy-formation model that 
we consider in this paper, and predicts that galactic spheroids and 
MBHs should grow earlier than in the "standard" model, which 
might result in larger event rates for LISA/SGO or eLISA/NGO (at 
least in the case of light MBH seeds). We will study the predic- 
tions of this alternative galaxy-formation model in detail in a future 
paper. 
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